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Abstract 

m 

Context. In 1992 we began a precision radial velocity survey for planets around solar-like stars with the Coude Echelle Spectrograph 
and the Long Camera (CES LC) at the 1.4 m telescope in La Silla (Chile) resulting in the discovery of the planet ( Hor b. We have 
continued the survey with the upgraded CES Very Long Camera (VLC) and the HARPS spectrographs, both at the 3.6 m telescope, 
until 2007. 

Aims. In this paper we present additional radial velocities for 31 stars of the original sample with higher precision. The observations 
cover a time span of up to 15 years and permit a search for Jupiter analogues. 

Methods. The survey was carried out with three different instruments/instrument configurations using the iodine absorption cell and 
the ThAr methods for wavelength calibration. We combine the data sets and perform a joint analysis for variability, trends, and peri- 
odicities. We compute Keplerian orbits for companions and detection limits in case of non-detections. Moreover, the HARPS radial 
velocities are analysed for correlations with activity indicators (Call H&K and cross-correlation function shape). 
Results. We achieve a long-term RV precision of 15m/s (CES+LC, 1992-1998), 9m/s (CES+VLC, 1999-2006), and 2.8 m/s 
(HARPS, 2003-2009, including archive data), respectively. This enables us to confirm the known planetary signals in ( Hor and 
HR 506 as well as the three known planets around HR 3259. A steady RV trend for e Ind A can be explained by a planetary com- 
panion and calls for direct imaging campaigns. On the other hand, we find previously reported trends to be smaller for /? Hyi and 
not present for a Men. The candidate planet e Eri b was not detected despite our better precision. Also the planet announced for 
HR 4523 cannot be confirmed. Long-term trends in several of our stars are compatible with known stellar companions. We provide a 
spectroscopic orbital solution for the binary HR 2400 and refined solutions for the planets around HR 506 and ( Hor. For some other 
stars the variations could be attributed to stellar activity, as e.g. the magnetic cycle in the case of HR 8323. 

Conclusions. The occurrence of two Jupiter-mass planets in our sample is in line with the estimate of 10% for the frequency of giant 
planets with periods smaller than 10 yr around solar-like stars. We have not detected a Jupiter analogue, while the detections limits 
for circular orbits indicate at 5 AU a sensitivity for minimum mass of at least lM Jllp (2Mj up ) for 13% (61%) of the stars. 



Key words, stars: general - stars: planetary systems - techniques: radial velocities 



1. Introduction dial velocity (RV) technique. Interestingly, many hot Jupiters 
jj] 1 have been found, a consequence related to the fact that the 
Cj ; The search for extra-solar planets has so far revealed approx- RV met hod as well as the transit method is more sensitive to 
imately 850 exoplanetfQ most of them discovered by the ra- short period p i ane ts. Out of 850 planets discovered so far, 65% 
have a period shorter than 1 year. Before the discovery of the 



* Based on observations collected at the European Southern first extra solar planet around a s olar-like star, the hot Jupiter 

Observatory, La Silla Chile, ESO programmes 50.7-0095, 51.7-0054, 51 Peg b (May or & Ouelozl lT995). it was widely expected that 

52.7-0002, 53.7-0064, 54.E-0424, 55.E-0361, 56.E-0490, 57.E-0142, planetary systems are in general similar to the solar-system and 

58.E-0134, 59.E-0597, 60.E-0386, 61.E-0589, 62.L-0490, 64.L-0568, this was also predi cted by most theoretical models as noted by 



66.C-0482, 67.C-0296, 69.C-0723, 70.C-0047, 71.C-0599, 072.C- iMarcv et al.l (120081). In the solar system. Jupiter is the dominant 

0513, 073.C-0784, 074.C-0012, 076.C-0878, 077.C-0530, 078.C-0833, planet amongst a n ot h er planets and causes the largest RV amp- 

ArcWve Facilif 80 ^ ^ ^ ° btamed ^ ^ SC16nCe litUde ' Therefore surve y s were set U P to search for P lanets with 

■ r , masses of 1 Mj UD and at distances of 5 AU from solar-like stars 

Radral velocrty data are avarlable in electronic form at the . i ,, r „ J V I N nnf t m, ■ j- T i •« 

™.„ (v . j . . - , 1 -,mn no cs ( e -g- Walker et al. 1995). The regime of Jupiter analogues is still 

CDS via anonymous ftp to cdsarc.u-strasbg.tr (130.79.128.5) or via v ° ; - " ° . ,\ . ° 

htt P ://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ sparsely explored because observations with long time-baselines 

1 http://exoplanet.eu 
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and precise RV measurements are required; e.g. Jupiter orbits the 
Sun in 12 years and induces an RV semi-amplitude of 12m/s. 

Th ere are many exopla n et search projects e.g. at Lick, 
AAT dO' Toole et alJ l2009t IWittenmver etail 1201 ll) Keck 
dCumming et alj 120081) . ELODIE/SOPHIE dNaef et alj [ 2005 ; 



Table 1. Targets with their sp ectral type SpT 
1991), visual magnitude V (Perrvman et al. 



Joffleit & Jaschekl 
19971) . distance d 



Bou chv et alj l2009al). CORALIE dSegransanetal.ll201oi) , and 



HARPS (Naefetal 
cision RV projects 



20101 iMavor et alJl201 lh . These high pre- 
have discovered a large fraction of the 
currently known planets and are continuously extending their 
time baselin es. Examples for discovered Jupiter analogues are 
GJ 777 Ab dNaef et ail Eool . a 1.33 A ft™, planet at 4.8A U3 
around a G6IV star, or HD 154345b d Wright et alj l2008t) . a 
0.95 Mj up planet at 4.5 AU around a G8V dwarf (all masses 
are Msini minimum masse s). Two more Jup iter-analogues 
were also recently reported bv lBoisse et al.l d2012l) : HD150706b 
(2.7M Jup , 7 AU) and HD222 155b (1.9 Mj up , 5.1 AU). 

The survey d escribed in this paper was begun in 1992 
dEndl et al.ll2002h with the Coude Echelle Spectrograph (CES) 
Long Camera (LC). With the advent of the CES Very Long 
Camera (VLC) in 1999, it was transferred to this instrument 
combination, and was continued in 2003 with the HARPS spec- 
trograph. The last observations for this programme were taken in 
September of 2007, although we have also made use of archival 
data acquired up to 2009. The survey covers a time span of up to 
15 years with RV precisions ranging fr om 15 m/s down to 2 m/s . 
A comparable survey was analysed by Witt enmver et alj {2006) 
and carried out in the northern hemisphere with the 2.7 m tele- 
scope at the McDonald Observatory. It started in 1988 with 24 
solar- like star j^] and 7 s ubgiants. When combined with CFHT 
data (Wal ker et alJI 1995b . it gave an even longer temporal cov- 
erage of up to 25 years, albeit with a somewhat lower precision 
(10-20m/s). 



2. The sample 

The o riginal sample of 37 solar-like stars was introduced in de- 
tail in En dl et al]d2002l) . Of these, the monitoring of six stars was 
stopped: HR 448, HR 753, HR 7373, Barnard's star, Proxima 
Centauri, and GJ 433. The first three had been observed tempo- 
rarily in 1996/97 as once promising metal rich targets, but were 
soon left out due to limited observing time. The latter three are 
M dwarfs which were included in a dedicated M dwarf survey 
with VLT+UVE S. For these stars re c ent an d more precise results 
are published in IZechmeister et all d2009l) . So we are left with 
the 31 stars listed in Table Q] along with some of their properties 
(spectral type, visual magnitude, distance, and stellar mass). 

All stars have a brightness of V < 6 mag and spectral types 
ranging from late F to K. There are two subgiant stars (J3 Hyi 
and 5 Eri) and two giant stars (HR 3677 and HR 8883) in the 
samplfl 

The sample includes six stars for which planet detec- 
tions have been claimed. These are HR 506 (Mayor et al.fl 



2 The planet was confirmed by Voet et al. (2005) who revised the 
semi-major axis to 3.9 AU and discovered a second inner planet (17. 1 d, 
0.057 Mj U p). 

3 There are three targets (S Eri, a For, and r Cet) common to both 
samples (we do not combine the measurement of both samples). 

4 HR 3677 a nd HR 8883 were indicate d in the Bright Star Catalogue 
as dwarf stars jHoffleit isT jaschek 1991). Therefore they entered our 
sample, however they are giants as indicated by their distances. 

5 HR 506b was announced by Mayor et al. at the 
XIX th IAP Colloquium in Paris (2003). We found 
no refereed publication. Information is available on 
http : //ob swww . uni g e . ch/ -udry/plane t/hd 10647. html 



(Ivan Leeuwenl 120071) . secular acceleration v r , and stellar mass 



M. 


Star 


alias 


SpT 


V [mag] 


d [pc] v r 


[m/s/yr] 


M[M Q ] 


HR77 


f Tuc 


F9V 


4.23 


8.59 


0.84 


1.06 [PM] 


HR98 


£Hyi 


G2IV 


2.82 


7.46 


0.86 


1.1 [D] 


HR209 


HR 209 


G1V 


5.80 


15.16 


0.01 


1.10 [G] 


HR370 


vPhe 


F8V 


4.97 


15.11 


0.16 


1.20 [G] 


HR506 


HR 506 


F9V 


5.52 


17.43 


0.02 


1.17 [G] 


HR509 


rCet 


G8V 


3.49 


3.65 


0.31 


0.78 [T] 


HR695 


k For 


GOV 


5.19 


21.96 


0.02 


1.12 [G] 


HR810 


iHor 


GOV 


5.40 


17.17 


0.06 


1.25 [V] 


HR963 


a For 


F8V 


3.80 


14.24 


0.17 


1.20 [G] 


HR1006 


Ret 


G2.5V 


5.53 


12.01 


0.61 


1.05 [G] 


HR1010 


<T 2 Ret 


G1V 


5.24 


12.03 


0.61 


1.10 [G] 


HR1084 


e Eri 


K2V 


3.72 


3.22 


0.07 


0.85 [DS] 


HR1136 


S Eri 


K0IV 


3.52 


9.04 


0.12 


1.23 [PM] 


HR2261 


a Men 


G6V 


5.08 


10.20 


0.01 


0.95 [G] 


HR2400 


HR 2400 


F8V 


5.58 


36.91 


0.02 


1.20 [G] 


HR2667 


HR 2667 


G3V 


5.56 


16.52 


0.06 


1.04 [G] 


HR3259 


HR 3259 


G7.5V 


5.95 


12.49 


0.30 


0.90 [G] 


HR3677 


HR 3677 


G0III 


5.85 


196.85 


0.00 


2.1 [G] 


HR4523 


HR 4523 


G3V 


4.89 


9.22 


0.53 


1.04 [G] 


HR4979 


HR 4979 


G3V 


4.85 


20.67 


0.07 


1.04 [G] 


HR5459 


a Cen A 


G2V 


-0.01 


1.25 


0.42 


1.10 [P] 


HR5460 


o-CenB 


K1V 


1.35 


1.32 


0.40 


0.93 [P] 


HR5568 


GJ 570 A 


K4V 


5.72 


5.84 


0.54 


0.71 [G] 


HR6416 


HR6416 


G8V 


5.47 


8.80 


0.22 


0.89 [G] 


HR6998 


HR 6998 


G4V 


5.85 


13.08 


0.01 


1.00 [G] 


HR7703 


HR 7703 


K3V 


5.32 


6.02 


0.37 


0.74 [G] 


HR7875 


4> 2 Pav 


F8V 


5.11 


24.66 


0.24 


1.1 [PM] 


HR8323 


HR 8323 


GOV 


5.57 


15.99 


0.04 


1.12 [G] 


HR8387 


e Ind A 


K4.5V 


4.69 


3.62 


1.84 


0.70 [G] 


HR8501 


HR 8501 


G3V 


5.36 


13.79 


0.19 


1.04 [G] 


HR8883 


HR 8883 


G4III 


5.65 


101.32 


0.00 


2.1 [G] 



References for mass esti mates : [D] iDravins et all ((1998), [DS] 
iDrake & Smith! d!993l). 1G1 iGravl <1988l) . [PM] Porto de Mello , priv. 
comm., [P] iPourbaix et alJ d2002h . Tri iTeixeira et all d2009h . [V] 
IVauclairetalJd2008l) . 

Table 2. The three used instruments/configurations with their 
wavelength reference, chosen spectral coverage and resolving 
power, and telescope diameter. 



Spectrograph 


Ref. 


A [A] 


R 


Tel. 


CES + LC 


h 


5367 -5412 


100000 


1.4m 


CES + VLC 


h 


5376-5412 


220000 


3.6 m 


HARPS 


ThAr 


3800 - 6900 


115 000 


3.6 m 



t Hor dKurster et al.ll2000h . e Eri (lHatzes et al.ll2000h . HR 3259 
(lLovis et al.ll2006l). HR 4523 dTinnevetal.1120111) , and recently 



a Cen B dDumusque et al In Sect. 5 we provide more de- 

tailed information on individual objects and we will stress those 
planet hypotheses. 



3. Instruments and data reduction 

We used three high resolution spectrographs that are briefly de- 
scribed below with more detail provided for the less known 
VLC+CES. Table|2]gives an overview of some basic properties 
of the three instruments. 
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3.1. CES + Long Camera 

In 1992 the survey started (1992-1 1- 03 to 1998-0 4-04) with the 
Coude Echelle Spectrograph (CES; lEnardl[l982l) and its Long 
Camera (LC) fed by the 1.4 m Coude Auxiliary Telescope (CAT) 
at La Silla (Chile). The CES+LC had a chosen wavelength cov- 
erage of 45 A and a resolution of 100000 (Tabled. A 2k x 2k 
CCD gathered part of one spectral Echelle order. An iodine gas 
absorption cell provided the wavelength calibration. More de- 
tails about the instrum ent, data analysis, as well as the obtained 
results can be found in Endl et a Dd2002|). Table[3]lists the radial 
velocity results. The median rms is 15.2 m/s when excluding the 
giants and targets with companions and trends as commented in 
Tableland reflects the typical precision. 

3.2. CES + Very Long Camera 

The Very Long Camera (VLC; iPiskunov etaD 1 19971) of the 
Coude Echelle Spectrograph (CES) was commissioned at the 
ESO 3.6 m telescope in La Silla (Chile) in April 1998 and de- 
commissioned in 2007. The VLC was an upgrade of the CES that 
doubled the resolving power to R — 220 000 - 235 000 as well as 
the CCD length so that 80% of the spectral coverage compared 
to the LC was retained (cf. Table[2]). This upgrade together with 
improved internal stability, and also the larger telescope aperture 
promised an improvement of the RV precision. For our sample 
we collected VLC spectra from 1999-1 1-21 to 2006-05-24. 

The VLC was fed by a fibre link from the Cassegrain fo- 
cus of the 3.6m telescope. A modified Bowen-Walraven image 
slicer provided an efficient light throughput at the high resolving 
power. It redistributed the light from the fibre with a 2" aperture 
via 14 slices to an effective slit width of 0.16" and resulted in a 
complex illumination profile in the spatial direction, i.e. perpen- 
dicular to the dispersion axis (Figs. Q] and 0. The right half of 
a 4 k x 2 k EEV CCD recorded part of one spectral order with 
the wavelength range of 5376-5412 A. In 2000-06-15 CCD#59 
was replaced by CCD#61 and in 2001-1 1-23 the CES fibre was 
exchanged. 

The CES + VLC employed the same iodine cell as the 
CES+LC for wavelength calibration. This cell was controlled 
at a temperature of 50°C. The RV modelling (Sect. \3A\ re- 
quires a high resolution and high signal-to-noise iodine spec- 
trum to reconstruct confidently the instrumental line profile (IP) 
of the spectrograph. In November 2008 we obtained a laboratory 
spectrum for our iodine cell with R = 925 000 and S/N-1000 
using a Bruker IFS125HR high-resolution Fourier Transform 
Spectrometer (FTS) at Lund Observatory. A filament lamp was 
used as a background source. To limit the light from adjacent 
wavelength regions a set of filters were applied: a coloured glass 
filter type VG1 1, a notch-filter to reduce the internal laser light, 
and a hot mirror to suppress the red spectrum. While the iod ine 
cell laboratory spectrum previously used by End l et al.l (120021) in 
their analysis of the CES LC data had only R - 400 000, the new 
scan ensures an iodine spectrum with a resolution almost 5 times 
higher than the resolution of the CES + VLC. 

The following properties of the CES + VLC spectra must be 
considered in the data analysis: The VLC spectra are contamina- 
ted by a grating ghost located in the middle of the CCD (Fig. [TJ 
and suffered also from stray light produced by the image slicer. 
Ripples are visible in the continuum of high S/N (~1000) spectra 
caused by interference in the fibre. This can be seen, for instance, 
in flatfield exposures. Also visible in flats are less efficient rows 
on the chip every 512 pixels, due to a smaller pixel size result- 
ing from the manufacturing process, which affect the wavelength 




spatial axis 



Figure 1. A small section of a VLC raw spectrum (star seen 
through the iodine cell; broad absorption features are stellar lines 
whereas narrow lines were introduced by the molecular iodine 
gas). The 14 slices span 400 pixels. The bright feature to the left 
near pixel row 2300 is a grating ghost. The area on the left side 
has a lower bias. The readout register is in the lower right corner 
of the chip (parallel clocking down, serial clocking to the right). 
Deeper stellar lines near row 2300 and 2700 have tails to the 
left and right. (The intensity scale is non-linear to bring out the 
discussed effects.) 



\ if 


k : 





200 400 600 800 1000 

spatial pixel 

Figure 2. Spatial profile of a VLC spectrum (cross-section for 
the rows 2400-2410 of the raw spectrum in Fig.[TJ linear intens- 
ity scale). 



solution. Moreover, as a peculiarity of the CES CCD electronics, 
a lower bias level is observed to the left of the spectra, caused 
by an electronic offset that occurs after processing a strong sig- 
nal. This effect is attributed to the video amplifier electronics 
and requires the readout of several CCD rows to properly dis- 
charge (P. Sinclair, ESO, 2011, priv. comm.). Hence subsequent 
CCD rows are affected which may cause systematic spectral line 
asymmetries and RV shifts depending on the spectral line depth. 
Moreover, since the iodine lines are weaker, they may not re- 
ceive the same shift as the stronger stellar lines and cannot cor- 
rect completely for this effectQ 

The spatial profile has a width spanning more than 400 pixels 
offering a large cross-section for cosmic ray hits (so-called cos- 
mics). For this reason the observing strategy aimed at three con- 
secutive spectra in one night to be able to identify cosmics as 
outliers. However, we did not use this cosmics detection method 
because cosmics could also be efficiently identified as deviations 
from the spatial profile in the optimum extraction. 



6 This effect looks similar to charge t ransfer inefficiency (CTI), which 
can also cause RV shifts of several m/s iBo uchv et alj20 09bl. However, 
CTI is caused by local defects on the CCD itself. 
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Table 3. Radial velocity results for all targets. For each instrument configuration the rms is calculated independently. RV data are 
not binned. Comments are on multiplicity (see also Table |5]for more information) and also indicate giants. 



Star 




CES + LC 






CES + VLC 






HARPS 




Comment 




N 


T 


rms 


ARV 


N 


T 


rms 


ARV 


N 


T 


rms 


ARV 








[d] 


[m/s] 


[m/s] 




[d] 


[m/s] 


[m/s] 




[d] 


[m/s] 


[m/s] 




<"Tuc 


51 


1888 


19.8 


14.5 


48 


2104 


9.2 


7.9 


1019 


2206 


1.9 


0.4 






157 


1887 


22.9 


18.6 


46 


1920 


7.7 


12.5 


2860 


1837 


2.6 


0.4 




HR 209 


35 


1572 


17.2 


17.6 


36 


1941 


11.1 


9.8 


48 


1401 


8.6 


0.5 




vPhe 


58 


1926 


15.6 


14.6 


35 


1910 


10.1 


8.2 


63 


1415 


2.7 


0.6 




HR506 


23 


1573 


28.0 


20.2 


42 


1910 


18.7 


17.7 


119 


1401 


10.7 


1.0 


planet 


rCet 


116 


1888 


11.4 


13.3 


61 


1920 


8.1 


8.8 


5373 


2244 


1.4 


0.3 




k For 


40 


1890 


722.9 


12.8 


45 


2094 


1134.0 


10.2 


74 


1401 


657.6 


0.5 


SB1 


l Hor 


95 


1976 


51.3 


16.9 


122 


2186 


37.1 


12.6 


1861 


1401 


13.1 


0.9 


planet 


a For 


65 


1889 


42.2 


31.6 


39 


1856 


16.7 


14.7 


191 


1401 


8.5 


0.7 


trend 


C Ret 


14 


184 


17.0 


14.2 


42 


1857 


15.3 


10.3 


63 


1401 


8.1 


0.5 




I 2 Ret 


58 


1976 


18.9 


14.6 


43 


1857 


10.1 


9.3 


82 


1411 


2.8 


0.5 




e Eri 


66 


1889 


12.2 


9.0 


69 


2186 


10.0 


8.1 


521 


1390 


5.5 


0.3 


planet(?) 


6 En 


48 


1888 


12.5 


11.7 


42 


1856 


7.5 


7.0 


152 


1410 


2.8 


0.2 




a Men 


46 


1852 


9.8 


10.7 


77 


2368 


8.4 


9.8 


188 


1308 


2.6 


0.3 




HR 2400 


53 


1924 


275.1 


23.0 


54 


2039 


523.8 


14.4 


63 


1296 


190.3 


0.8 


SB1 


HR 2667 


66 


1934 


15.1 


18.5 


64 


2329 


7.5 


11.4 


68 


1296 


1.6 


0.4 




HR 3259 


35 


1851 


16.5 


11.3 


61 


2367 


9.2 


7.7 


435 


1294 


3.7 


0.3 


three planets 


HR 3677 


34 


1924 


492.4 


15.1 


38 


2044 


1253.9 


8.2 


66 


1287 


870.4 


0.5 


SB1, giant 


HR 4523 


27 


1925 


14.9 


12.2 


57 


2276 


6.7 


8.9 


253 


1608 


3.4 


0.3 




HR 4979 


52 


1933 


11.7 


10.8 


58 


2329 


9.2 


10.2 


460 


1286 


3.7 


0.3 




a Cen A 


205 


1852 


166.5 


10.7 


1074 


1929 


97.7 


10.2 


5029 


1206 


21.0 


0.2 


SB1 


aCenB 


291 


1852 


203.5 


9.3 


54 


1770 


247.6 


7.7 


255 


1206 


191.9 


0.2 


SB1 


GJ 570 A 


40 


384 


6.9 


11.4 


87 


2284 


10.2 


6.5 


47 


1853 


2.7 


0.3 




HR6416 


57 


1845 


23.8 


12.6 


59 


2278 


23.4 


9.4 


76 


1310 


7.5 


0.4 


trend 


HR 6998 


51 


1789 


15.3 


20.8 


23 


2062 


9.6 


9.8 


70 


1044 


1.2 


0.4 




HR 7703 


30 


1042 


10.3 


11.6 


31 


2039 


7.6 


8.0 


79 


1735 


4.8 


0.4 


trend 


<fr Pav 


90 


1969 


32.1 


25.6 


200 


2062 


17.1 


23.9 


63 


2158 


4.0 


1.2 




HR 8323 


20 


1067 


14.4 


14.6 


31 


2124 


11.4 


10.1 


318 


1413 


3.7 


0.5 




e Ind A 


73 


1888 


11.9 


9.1 


54 


2124 


7.1 


7.9 


457 
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58 


1413 


15.2 


0.6 


trend 


HR 8883 


31 


1258 


63.1 


31.3 


30 


2125 


66.9 


23.6 


45 


1401 


61.7 


1.6 


giant 



Listed are the number of observations A', the time baseline T, the weighted rms of the time series and the effective mean internal radial velocity 
error ARV. 



The VLC spectra were reduced with standard IRAF-tasks in- 
cluding subtraction of the overscan and a nightly master-bias, 
2D flat-fielding, scattered light subtraction, and optimum extrac- 
tion (lHornelll986l) which also removes cosmics. The scattered 
light was defined left and right of the aperture with a low- 
degree polynomial used to interpolate across the aperture. This 
was done row-by-row and afterwards "smoothed" in the disper- 
sion direction with a high-order spline to account for the above 
mentioned features and then as scattered light subtracted from 
the spectra. Finally, the science spectra were roughly calibrated 
with a nightly ThAr spectrum to provide an initial guess for the 
wavelength solution which is later refined with the iodine spec- 
trum in the subsequent modelling process. The whole data re- 
duction process largely removed the artefacts described above, 
however residual deviations are likely to still exist in the RVs of 
the VLC data. The typical precision is 9.4 m/s calculated as the 
median rms in Table|3]for the stars without comments. 

3.3. HARPS 

With HARPS we monitored our targets from 2003-11-06 to 
2007-09-21 (2009-12-19 including archive data). The HARPS 
spect r ograph is desc ribed in the literature (e.g. iMavor et al.l 
2003; Perje et al. 2004). It is fibre fed from the Cassegrain focus 
of the 3.6 m telescope and located in a pressure and temperature 
stabilised environment. An optical fibre sends light from a ThAr 



lamp to the Cassegrain adapter for wavelength calibration. For 
the RV computation via cross-crorelation with a binay mask 72 
Echelle orders ranging from 3800 A to 6900 A are available, a 
region much larger than for the CES. 

We made use of the ESO advanced data products (ADP) to 
complement our time series which sometimes also extended the 
timebase. This archive provides fully reduced HARPS spectra 
including the final radial velocities processed by the pipeline 
DRS03.5 (data reduction software). The radial velocities are cor- 
rected for the wavelength drift of the spectrograph (if measured 
by the simultaneous calibration fibre) and the RV uncertainty 
estimated assuming photon nois^J The mean RV uncertainties 
range from 0.2 to 0.8 m/s and do not include calibration errors, 
guiding errors, and residual instrumental errors. For data ana- 
lysis a stellar jitter term (> 1.6 m/s) will be added in quadrature 
(see Sect.©. 

We recomputed with the HARPS DRS some of these 
archival RVs that suffered in the cross-correlation process from 
a misadjusted initial RV guess (off by more than 2 km/s) or from 
an inappropriate binary correlation mask. A different mask, e.g. 
K5 instead of G2, can produce RV shifts up to 20 m/s. The pub- 
licly available archive data originate from other programs such 



7 http : //www . eso . org/ sci/ f acilities/lasilla/ instruments/harp 

8 The pertinent information can be found in the *CCF_A.fits-file 



header (keywords RVC and DVRMS). 
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Figure 3. A VLC spectrum of t Cet (without iodine cell) and 
the arrangement of the 200 pixel chunks with their individual Q- 
factors (blue squares). The intensity maximum of the spectrum 
is set to unity. The intensity declines to the edges due to instru- 
mental effects (Blaze function). 



as short-term asteroseismology campaigns or the HARPS GTO 
(guaranteed time observations). The latter complemented our 
data with additional measurements, and in some cases provided 
a data set that outnumbers our own in terms of number of meas- 
urements and time base. We use only data taken in HARPS high 
accuracy mode (HAM), while we leave out data taken with iod- 
ine absorption cell or in high efficiency mode (EGGS, "Extra 
Good General Spectroscopy") which uses a different fibre, a dif- 
ferent injection method, and no scrambler and has a lower stabil- 
ity and a different zero point. Furthermore, spectra with a signal 
to noise of S/N < 50 are also left out. 

The HARPS data provide an absolute RV scale which is 
shown in Figs. [1914231 and serves as our references frame into 
which the other instruments are transferred. Note that relative 
RV measurements are more precise than the absolute RVs (i.e. 
more precise than accurate fj The median rms in Table [3] is 
2.8 m/s for stars without comments. Note that several of our stars 
are active, so that this value is higher than the precision of 1 m/s 
usually quoted for HARPS. 

To improve the combination of the HARPS and VLC data, 
spectra were taken in a few nights with both spectrographs im- 
mediately after each other making use of an easy switch possible 
with the common fibre adapter installed in May 2004. 

3.4. Details of the RV computation for the CES+VLC data 

To compute the R Vs of the VLC spe ctra we used the AUSTRAL 
code described in lEndl et al.l (2000) which is bas ed on the mo- 
delling technique outlined in iButler et all dl996l) . The spectral 
order was divided into 19 spectral segments (chunks) with a size 
of 200 pixels ( 1 . 8 A) which we empirically found to yield the op- 
timal RV precision. The reasons could be that for smaller chunks 
the stellar RV information content becomes too small. For a lar- 
ger chunk size, on the other hand, probably the assumption that 
the instrumental profile (IP) is constant over the chunk breaks 
down, or the discontinuities of the wavelength solution by the 



9 Systematics in absolute (spectroscopic) RVs arise e.g. from the stel- 
lar mask, gravitation al red shift and convective blue shift of the star (e.g. 
iPourbaix et alj2002h . 



mentioned smaller and less efficient pixel rows are more prob- 
lematic. 

The stellar spectrum can be shifted across the CCD by sev- 
eral pixels due to the barycentric velocity of the Earth (c alcu- 
lated with the JPL ephemerides DE200, e.g. IStandishll 1 9901) and 
offsets in the instrument setup. To ensure that the same stellar 
lines fell in the same chunk, we shifted the chunks to the proper 
spectral location according to the barycentric correction (1 pixel 
is ~ 500 m/s). So instead of having fixed chunk positions with 
respect to the CCD as originally implemented in the AUSTRAL 
code, this modification ensures always the same weighting factor 
for each chunk. The final wavelength solution in each chunk is 
provided by the iodine lines which record the instrumental drifts 
and offsets. 

Figure [3]illustrates the alignment of the chunks with respect 
to the stellar spectrum. This placement of the chunks tries to 
avoid splitting up stellar lines between adjacent chunks. As one 
can see, the chunks contain only a few deep stellar lines or some- 
times none. To quantify this, we calculated the quality factor Q 
(IConneslll985t IButler et aL | fl996t iBouchvet all 1200 ll) for each 
chunk in a stellar template! This factor sums in a flux- weighted 
way the squared gradients in a spectrurrFT, hence measuring its 
RV information content. For photon noise, the estimated RV un- 
certainty is inversely proportional to Q, i.e. ARV ~ ^. Hence, we 

weight each chunk RV with Q 2 when computing the RV mean. 
Chunks with Q < 7 000 were discarded (cf. Fig. [3] right axis). 
For comparison, the quality factor is 2=12857 for the whole 
spectral range in Fig. [3] and Q - 61 000 for an iodine spectrum 
(e.g. the spectrum of a featureless B-star taken through the iod- 
ine cell). 

For t Cet (GJ 71, HIP 8102, HR 509, HD 10700) which is 
known as an RV constant star (t he HARPS data h ave an rms 
scatter of -1.1 m/s HARPS data. iPepe etai1l201 ll this work), 
we achieve with the CES+VLC a long-term precision of 8. 1 m/s 
(Fig. [19] Table[3]l. The internal RV errors of the individual spec- 
tra (~8.8m/s), calculated as the errors of the mean RV of the 
chunks (rms/ V-VchunkX are of the same order as the rms of the 
time series implying a fair error estimation. 



3.5. Combining the LC and VLC data 

The problem of instrumental offsets, i.e. different radial velo- 
city zero points, occurs when data sets originate from different 
instruments (e.g. IWittenmver et al. 2006) or after instrumental 
changes/u pgrades. For instan ce, an offset of -1.8 m/s was re- 
ported by Rivera et alj (|2010) after upgrading the Keck/HIRES 
spectrograph with a new CCD . An offset of only 0.9 m/s was 
mentioned bv lVogt et al.l d2010h when combining Keck and AAT 
data. 

As described above we have used three different instru- 
ments/instrument configurations and we are also faced with the 
problem of the instrumental offset. There are basically two dif- 
ferent methods for combining the data sets: (1) Simply fitting 
the offset, i.e. the data sets are considered to be completely in- 
dependent and the zero points are free parameters in the model 
fitting. (2) If possible, measuring the offset physically by making 



10 This template was used in the modelling and obtained via deconvo- 
lution from a ste llar spectrum taken without iodine cell as described in 



Endl et alJ 



2002). 



ZA«)(Hf 



ZA(0 



where A(i) is the flux in the f'-th pixel. 



5 



Zechmeister et al.: The planet search programme at the ESO CES and HARPS 



use of some known relation between the data sets/instruments to 
keep the offset fixed. 

In fact, we can measure the offset for the LC and VLC data 
albeit with a limited precision. The LC and VLC spectra were 
taken through th e same iodine cell, i.e. the same wavelength cal- 
ibrator. Because lEndl et al.1 d2002l) calculated the LC RVs with 
different stellar templates and an iodine spectrum of lower res- 
olution than used in this work, we re-calculated the RVs for all 
LC spectra with the same VLC stellar template (which is shorter 
than the LC spectra) and the new iodine cell scan to have the 
same reference for the LC and VLC. The re-calculated RVs are 
verified to have a precision similar to the published LC data. 

Then we computed the mean of the re-calculated LC and 
VLC time series. If a star has a constant RV, one would expect 
that the means of both time series are the same, i.e. the offset 
RVvlc _ RVlc = within the uncertainties of the means (cr L c 
and o~ vlc)- This can be tested with the f-statistics, in particular 
Welch's f-test (for two independent samples with unequal sizes 
and variances). We suggest that keeping the offset fixed is valid, 
if the quantity 



t = 



RVlc - RVvlc 



with 



LC 



VLC 



N LC N vlc 



(1) 



is not rejected by the Null-hypothesis. The parameter s is an es- 
timate for the standard error of the difference in the means and 
is calculated from sample variances aj and sample sizes N,. The 
variable t follows a f-distribution with v degrees of freedom^. 
For instance, for \ t\ > 1.7 and v > 30 the difference in the means 
is significant with a false alarm probability of FAP<10%. For 
some of our stars the FAP for the offset difference is not signi- 
ficant: 5 Eri (64%), e Eri (23%), HR 209 (92%), and HR 3259 
(15%). However, from Fig. [4] it can be seen that there are also 
stars having significant offsets leaving doubts whether the offset 
can be kept fixed in general. 

Figure [4] shows that for our sample an average offset of 
10.4 m/s (+7.7 m/s) remains for the RV constant stars when com- 
paring the RV means of the VLC and LC data. This offset might 
be due to systematics in the deconvolution process of the stel- 
lar template or in the modelling. For example, due to the dif- 
ferent resolution, the LC data have to be modelled with a dif- 
ferent chunk size (154 pixels to cover two VLC chunks). We 
corrected all re-calculated LC RVs for this systematic offset. 
Finall y, we adjusted th e RV mean of the published LC time 
series dEndl et a l. 2002) to fit the R V mean of the re -calculated 
time series. In Figs. [19H23] the LC dEndl et al.ll2Q02h and VLC 
data are always shown relative to each other with the measured 
and corrected offset (and not with a fitted RV offsets that could 
have been taken from our fit results presented below in Sect. |4]l 
to conserve the true measurements. 

The uncertainty of the offset found in the sample is rather 
large for it to be considered a fixed value. On the other hand 
the approximately known offset can hold important information, 
in particular in the case of HR 2400 or e Ind A. Therefore we 
choose a compromise between a fixed offset and a free offset 
when fitting a function. Because one expects the difference of the 
zero point parameters to be zero (cvlc - <?lc ~ 0), we introduce 
in the ^-fitting a counteracting potential term rf (also called 



The 



effective 



(st/Nt + s\lN 2 f 



degree 



tfjNtfm - 1) + (s 2 2 /N 2 y/(N 2 - 1) 

the sample variances. 



of 
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Figure 4. Difference between the means of the VLC and re- 
computed LC time series for all our stars ordered by spectral 
type. For RV constant stars (black filled circles) there occurs a 
systematic offset of 10.4±7.7m/s (red solid line and red dashed 
lines). Stars shown with open circles were not included in the 
offset analysis (for reasons see comments in TableO. The shown 
error bars correspond to the uncertainty in the means, i.e. para- 
meter s from Eq. ([TJ. 



penalty function, e.g. Shpor er et alJl2010h . that increases when 
the zero point difference becomes larger 



X 



2 , 2 

■X +1 



with 



CVLC - C L C 



(2) 



The resulting x 1 (when minimising x 1 ) will be higher com- 
pared to that obtained when fitting with free offsets but lower 
than for fixed offsets. The parameter s determines the coupling 
between the offsets. After performing the fit it can be checked, 
if the fit has spread the zero points too much (if r\ » 1 or if 
there are large jumps in the model curves in Figs. |T9l423l >. For s 
we attributed the uncertainty of the offset correction of 7.7 m/s 
leading to a weak coupling. 

It is worth mentioning, that in Bayesian analysis x 1 can be 
identified with the likelihood when assuming a Gaussian dis- 
tribution for the prior information that the expected zero point 
difference is zero. 



3.6. Combining the CES and HARPS data 

In principle, VLC and HARPS data could be combined in a simi- 
lar way. They have different wavelength calibrators, but, since 
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there are some nights with almost simultaneous observations 
(within minutes), they are closely related in time. The difference 
between these consecutive measurements should be zero so that 
it is tempting to bind directly the time series by means of those 
nights. However, this does not account for fluctuations due to 
the individual uncertainties. Again a coupling terrrP^I would be a 
more secure approach. 

However, for reasons of simplicity we choose a fully free 
offset between the HARPS and the CES data. Because the VLC 
and the HARPS time series overlap well this is less critical, in 
contrast to the LC and VLC time series which are separated by 
a 2-year gap. The relative offsets between the CES and HARPS 
data as illustrated in Figs.[T9l-l23lcorrespond to the common best 
fitting model (constant, slope, sinusoid, or Keplerian; cf. Sect. [4]). 

4. Analysis of the radial velocities 

In this section we describe our data analysis and the general re- 
sults of the survey, while some individual objects are discussed 
in detail in Sect. [5] The tests which we perform hereafter were 
repeated on the residuals of the binaries and planet hosting stars 
to search for additional companions and are indicated as objects 
with index r in Table [4] 

4.1. Preparation of RV data and jitter consideration 

Before the data analysis we binned the data into 2-hr intervals by 
calculating weighted means for the temporal midpoint, RV, and 
RV error. The 2-hr interval will especially down-weight nights 
from asteroseismol ogy campaign s (see Sect. 13 - 3 b and reduce the 
impact of red noise (Baluev 2012), while resulting in nightly av- 
erages for most other nights and still permitting to search for 
planets with periods as short as o ne day. Such interv als are also 
employed for solar-like stars (e.g. lRivera et alJuOlOl) to average 
out the stellar jitter, i.e. intrinsic stellar RV variation caused by, 
e.g., oscillation or granulation in the atmospheres of the stars. 
While the Sun has an oscillation ti mescale of ~5 min, its gran- 
ulation^ has lifetime <25min (see iDumusque et all 1201 lbl for 
adequate observing strategies). However note, that in our own 
survey we have usually taken three consecutive spectra in one 
night covering in total only 5-10 min, which is not sufficient to 
average out all those intrinsic stellar RV variations. To investig- 
ate the short-period jitter, we calculated the weighted scatteiP^I in 
each 2-hr bin with at least 2 measurements and then the weighted 
mean of these scatter^. Table [A~2l lists the jitter estimate <7jj t;T 
from the HARPS data for each star and the mean time scale t 
accessible for this estimate within the 2 h bins. Note that these 
time scales may not sufficiently cover the real jitter time scale in 



13 For one simultaneous measurement taken at the time t„ this term 

ch - Cvlc , . 

— where c is the zero 



could be written as i] = — - 

^ARVn(t s y-+ARV VLC (t s ) 2 
point parameter, ARV(t s ) are the individual errors of the simultaneous 
measurements, and the indices H (HARPS) and VLC indicate the in- 
struments. The VLC time series must be a priori adjusted by a zero 
point such that RV yLC (t s ) := RV H (t s ). 

14 There is also meso- and supergranulation (life times up to ~ 
30 hour s) which take place on different size scales Dumusque et al. 
d2011bh . 

15 Weighting of the i-th measurement with its internal error w t ~ 

16 Weighting of the j-th bin with the number of measurements nj and 
the mean internal error in that bin: w ; ~ n -Ja 1 , .. Note that bins with 

J J 1 lnt,y 

more measurements usually cover larger time intervals and get more 
weight. 



all cases. Therefore these estimated jitter values were not used 
in a further analysis. 

There can be also a long-term jitter with time scales of few 
days to weeks related to the rotation period (due to the appear- 
ance and disappearance of spots) or up to few ye ars due to the 
magnetic cycle of a star. llsaacson & Fischer! ( 20 10) provide jitter 
estimates as a function of B — V colour and chromospheric activ- 
ity index Shk based on Keck observations for more than 2600 
main sequence stars and subgiants. Using these relations and the 
median Shk values in the HARPS data (which in most cases 
agree well with other literature values; see Table IA.2I ). we es- 
timate the jitter crj it i ong for our stars (Table lAT2l . The jitter terms 
are usally > 2m/s. For GJ 570 A and e Ind A the expected jit- 
ter is only 1.6 m/s. Both are K dwarfs with 1 .0 < B - V < 1.3 
and, according to Isaac son & Fischeil d2010l) . those stars have 
the lowest level of velocity jitter decoupled from their chromo- 
spheric activity. The jitter terms <Tj it long were added in quadrat- 
ure to the internal errors for all stars and lead to a more balanced 
fit with the CES data. Morevover, to cross-check whether detec- 
ted RV signals might be caused by those kinds of stellar activity 
we will also analyse in the HARPS data correlations between 
the RV data and activity indicators such as Ca II H&K emis- 
sion and variations of the bisector (BIS) and the FWHM of the 
cross-correlation profile (Sect. l4.5t . All RVs and HARPS activ- 
ity indicators are online available. 

In fitting the data we accounted for the secular acceleration 
of the RVs (as given in Table [TJ. This perspective effect can be- 
come a measurable positive trend in some high proper motion 
stars dSchlesingerj|1917t iKurster et ail 120031: IZechmeister et all 
2009). In our sample, e Ind A has the highest secular accelera- 
tion with 1.8m/s/yr. Its contribution is depicted in Fig.|23]by a 
dashed line. 



4.2. Excess variability 

To investigate objects for excess variability it is common to com- 
pare the observed scatter with a noise estimate. A significantly 
larger scatter indicates variability. Because internal errors ARV, 
and jitter estimations <7ji t are available, the quality of each mea- 
surement is assessed and allows us to weight the measurements 
in the ^-statistics, w,- = -K = „ „„' , _ 2 ■ As the scatter we calcu- 



ARVUcr: 

late the weighted rms which is here defined as 



N 1 



Z 



N 



X 1 



N-vW 



(3) 



where W = Y* w i is the sum of the weights and v the num- 
ber of model parameters. Outliers with a large uncertainty will 
contribute less to the rms. The factor t^- is a correction that 

N—v 

converts the uncorrected and biased variance into an unbiased 
variance, i.e. to have an unbiased estimator for the population 
varianc^H In the unweighted case (w, = 1, W = N) we 
obtain the well known formula for the unbiased rms: rms = 



17 Note however, that the square root of this variance, rms = 
\/rmc2 i s no t a unbiased estimate of the population standard deviation 



dDeakin & Kildealll999h . 
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Furthermore, we define the weighted mean noise term cr via 
the mean of the weight^ 



cr = 




(4) 



Again lower-quality measurements will contribute less to the 
mean noise level. 

With these definitions the reduced^ 2 can be easily expressed 
as the ratio of weighted rms to weighted mean noise level 



2 2 

2 x rms 

^red 



N- v 



(T- 



(5) 



To test for excess variability we have to fit a constant and to 
calculate the scatter around the fit. For the joint analysis we ac- 
count for the zero point parameter of each data set when fitting 
a constant as outlined in Sect. 13.51 Note that the probability for 
the excess variability Prob(^- 2 d ) is directly reliant on a proper 
estimate for the noise level cr. Also note that the tests in the 
next sections employ model comparisons and the jitter estim- 
ate enters only indirectly through fitting with modified weights. 
Table|4]summarises for the combined data set the weighted noise 
term cr, the weighted rms, and the x 2 -probability for this test. 
Table IaTI lists additionally the individual rms (columns 5-7 la- 
belled rms constant ) for each instrument. These values can differ 
from Table[3] because in Table |4]secular acceleration is accoun- 
ted for, jitter has been added, the data are binned, and the LC and 
VLC offsets are coupled. Because the HARPS data have a much 
higher precision, they dominate the statistics. 

The small ^-probabilities for most of the stars indicate that 
they are variable with respect to our noise estimate cr. However, 
10 stars (and the RV orbit residuals of 8 stars) have Prob(y 2 ) > 
1%, i.e. they show only low or no excess variability. In five cases 
the scatter is smaller than the noise level, i.e. ^- 2 d < 1, imply- 
ing an overestimation of the noise level. Indeed, for four stars 
(£ Tuc, t Cet, HR 2667, and HR 6998) the jitter estimate cr jiUong 
in Table \A2\ is higher than the scatter of the HARPS measure- 
ments (Table lATt . The reason for jitter overestimation might be 
a somewhat lower precisi on of the Keck sample from which the 
jitter relation was derived dlsaacson & Fischeii20Toh . 

4.3. Long-term trends 

Because potential planets or companions can have orbital peri- 
ods much longer than our observations, these objects may betray 
themselves by a trend in the RVs. We searched for trends by fit- 
ting a slope to the data and derived its significance via the fit 
improvement with respect to the constant model (previous Sect.) 
via 

L slope 



f slope =(JV" 4)- 



2 _ 2 

. ^constant X$\ 



x. 



slope 



or when expressed with unbiased weighted variances 



(Af-3)rms 2 -(JV-4)rms 2 



"* slope 



slope 



(6) 



(7) 



slope 



The associated probability for this F-value follows a F\ji-\- 
distribution (4 parameters: 1 slope, 3 zero points). Again Table|4] 



18 Another point of view leads to the same result: Gaussian errors are 
added in quadrature. Hence the trivial weighted mean is J-^ 2 = 



summarises the test for long-term trends. When adopting a false 
alarm probability threshold of < 1CT 3 fitting a slope improves 
significantly the rms of all binaries as well as that of p Hyi, 
t Cet, GJ 570 A, e Ind A and the residuals of HR 506. We note 
that for p Hyi, a For, GJ 570 A, HR 6416, HR 7703, e Ind A, 
and HR 8501 the trend is a sufficient model (regarding sinusoid 
and Keplerian fit, see next Section), because of the smaller FAP 
or weighted rms (i.e. smaller ;f 2 ed ). For these stars the trend is 
depicted in Figs. [T914231 

Some of our stars have known wide visual companions with 
a known separation p listed in Table[5] Whether these objects are 
able to cause the observed trend, can be verified by the estimate 

cm 



m n smi > > = 5.6 • 10 M-, 



\ZA\ 



J up 



m/s/yr VAU7 



with the radial acceleration |zaI ^ I^aI sin i of the observed com- 
ponent A and the projected separation p < r between both 
components. For comparison, Jupiter at 5.2 AU can accelerate 
the Sun by 6.6 m/s/yr. Table [5] summaries the information about 
known wide companions and shows that the minimum compan- 
ion masses derived from the measured slopes are below 0.5 M Q 
for a For, GJ 570 A, HR 6416, HR 7703, and HR 8501. These 
masses are in agreement with the masses as expected from the 
spectral type of their companions. However, e Ind B cannot ex- 
plain the trend seen for e Ind A (see Sect.[5]for details). 

The other possibility for a trend is an unknown and unseen 
companion. Whether the strength and the long duration of a 
trend is still compatible with a planetary companion, can be es- 
timated more conveniently, when Eq. © is expressed in terms 
of the orbital period P which is also unknown but has to be (for 
circular orbits) at least twice as large as the time span of obser- 
vations T. With Kepler's 3 rd law £ = Gf?, Eq. © can be 



written as 



\2/3 



/ P \ 4/3 

m B sin / > G~ 1/3 |z| - (M A + m B ) 2 

m/s/yr \yrj \ M / 

For e Ind A we find its companion to have m sin i > 0.97 Mj up 
forP>30yr. 

4.4. Search for Periodicities and Keplerian orbits 

To search for the best-fitting sinusoidal and Keplerian orbits, we 
employed the generalised Lomb-Scarg le (GLS) algorithm de- 
scribed in lZechmeister & Kiirsterl d2009l) . It was adapted to treat 
all three data sets with different offsets and also incorporates the 
weak offset coupling described before (Sect. [3~5l l. Searching for 
sine waves is a robust method to find periodicities and orbits with 
low eccentricities, while for highly eccentric orbits the Keplerian 
model should be applied. 

Figures [5] Q2] [13] [14] and [15] show the periodograms for 
some of the stars discussed here. The periodograms are norm- 
alised as 



2 _ 2 

^constant ^sin 1 

P = 2 an " ^ Ke P 

-^constant 



2 _ 2 

/^constant ^Ke: 

y 2 

A constant 



(10) 



involving the x 2 of the constant, sinusoidal, and Kepl erian 
model, respectively. Analogo us to ICumming et al.l 120081 and 
IZechmeister & Kiirsterl d2009l) . we calculated the probabilities 
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Table 4. Summary of the tests for excess variability, slope, and periodicities (sinusoidal and Keplerian) for the combined data set. 
Significant trends and sinusoidal periods are printed in bold face (FAP< \%6). An index "r" to the name of the star (first column) 
indicates tests on the residuals as derived from the most significant model to the original data. 



Star 


N hia 


T 


<x 


rms 


S 




rms slopc 


slope 


FAP 


rms si „ 


p 

* sin 


FAP 


rms Kc p 


^Kep 


FAP 






[yr] 


[m/s] 


[m/s] 






[m/s] 


[m/s/yr] 




[m/s] 


[d] 




[m/s] 


[d] 




i Tuc 


197 


17.0 


3.08 


2.24 




1 


2.24 


0.04 


0.57 


2.18 


14.8 


1 


1.99 


551. 


3- 10-" 


ySHyi 


109 


16.0 


3.34 


3.84 




0.014 


3.29 


1.79 


8.9 10 9 


3.07 


4340. 


1.8 10 7 


2.96 


1.03 


5 10 8 


HR 209 


53 


14.0 


6.57 


9.44 


A 


1.4- 10~ 5 


9.28 


1.81 


0.1 


7.26 


1.24 


0.0064 


5.92 


2.52 


5-10- 6 


vPhe 


71 


14.9 


4.72 


5.23 




0.1 


5.17 


-0.96 


0.13 


4.59 


1.98 


0.3 


4.37 


1.09 


0.037 


HR 506 


60 


14.0 


4.68 


12.2 


P 





12.3 


-0.69 


0.73 


6.69 


963. 


8.7- 10 12 


6.44 


963. 


3-10-" 


HR 506 r 








6.53 


A 


2.5- 10~ 5 


5.89 


-2.41 


0.00041 


5.31 


3.60 


0.021 


4.68 


3.61 


0.00023 


r Cet 


339 


17.1 


2.41 


1.37 




1 


1.33 


-0.19 


5.2- 10 6 


1.28 


380. 


5.5 10 7 


1.25 


382. 


7-10^ 


k For 


78 


14.8 


3.69 


717. 


B 





49.4 


-701. 





8.13 


>30yr 


210 139 


4.03 


10700. 


5-10- 156 


k For r 








3.87 




0.18 


3.88 


0.29 


0.38 


3.50 


4010. 


0.67 


3.11 


3720. 


0.0018 


t Hor 


154 


14.8 


5.79 


31.7 


P 





31.8 


-0.11 


0.92 


15.4 


307. 


2.6 10 44 


12.9 


307. 


2-10~ 53 


( Hor, 








12.5 


A 





12.3 


-2.14 


0.029 


10.9 


5.72 


5.6 10 6 


9.83 


5.81 


6-10- 13 


a For 


75 


14.8 


4.28 


13.2 


B 





4.25 


-11.5 


7.710" 37 


4.25 


>30yr 


7.2- 10 32 


4.24 


8230. 


7-10- 30 


a For,. 








4.19 




0.58 


4.22 


0.03 


0.93 


3.77 


350. 


0.73 


3.62 


72.6 


0.37 


f 1 Ret 


48 


9.9 


5.09 


12.4 


A 


MO- 33 


11.5 


-4.14 


0.0059 


9.01 


21.5 


0.0013 


7.51 


21.4 


9-10- 6 


<T 2 Ret 


89 


14.8 


4.69 


5.23 




0.064 


5.13 


0.85 


0.037 


4.95 


2.48 


1 


4.41 


3.90 


0.01 


e Eri 


79 


14.8 


5.41 


8.22 


A 


8.1 • 10 10 


8.25 


0.54 


0.49 


6.85 


3.11 


0.0029 


6.02 


1.44 


4-10~ 6 


(5Eri 


75 


14.8 


2.78 


3.23 




0.027 


3.24 


-0.09 


0.61 


2.85 


3.07 


0.29 


2.69 


2.66 


0.045 


a Men 


102 


14.2 


3.23 


3.12 




0.66 


3.12 


-0.26 


0.27 


2.75 


880. 


0.0095 


2.61 


837. 


0.00051 


HR 2400 


77 


14.4 


4.53 


266. 


B 





215. 


-103. 


1.4 10 8 


107. 


5860. 


6.1 10 26 


5.89 


9490. 


3-10-" 2 


HR 2400, 








5.56 




0.0044 


5.38 


-1.46 


0.016 


5.23 


1.25 


1 


4.48 


1.88 


0.0028 


HR 2667 


88 


14.4 


4.82 


3.09 




1 


3.08 


-0.31 


0.26 


2.91 


3.49 


1 


2.63 


1.63 


0.03 


HR 3259 


191 


14.2 


2.49 


4.12 


P 


4.4- 10~ 32 


4.10 


0.58 


0.086 


3.21 


8.67 


1.210 17 


3.22 


8.67 


4-10-' 6 


HR 3677 


61 


14.4 


3.97 


1000. 


B 





257. 


808. 


1.6 10 35 


35.2 


>30yr 


7.8 lO 79 


8.68 


>30yr 


M0- 107 


HR 3677, 








8.39 




6.6- 10 27 


8.46 


0.16 


0.82 


6.16 


8.79 


6.2 10 5 


5.84 


31.2 


4-10~ 5 


HR 4523 


100 


15.3 


3.35 


3.47 




0.29 


3.49 


-0.01 


0.97 


3.32 


2.92 


1 


2.91 


1.04 


0.00082 


HR 4979 


151 


14.6 


2.78 


3.96 


A 


2- 10 12 


3.91 


-0.75 


0.026 


3.45 


1.00 


3.710 6 


3.25 


1.00 


6-10- 9 


a Cen A 


121 


14.2 


3.31 


101. 


B 





10.2 


131. 





3.03 


>30yr 


5.1 10 174 


2.84 


>30yr 


8-10- 173 


a Cen A,. 








2.79 




0.99 


2.81 


0.00 


1 


2.51 


89.1 


0.0066 


2.29 


179. 


2-10~ 6 


o-CenB 


82 


14.2 


3.96 


179. 


B 





19.1 


-154. 





5.06 


>30yr 


1.2 10 116 


4.92 


9370. 


MO- 113 


a Cen B, 








4.79 




0.055 


4.81 


0.18 


0.61 


4.08 


390. 


0.0091 


3.51 


387. 


1- 10 6 


GJ 570 A 


64 


12.0 


2.89 


5.12 


B 


2-10~ 15 


3.71 


-2.72 


3.5- 10 10 


3.74 


>30yr 


1.6 10 5 


3.46 


3390. 


2-10~ 6 


HR 6416 


77 


14.5 


3.99 


11.0 


B 





3.83 


9.40 


1.7 10 35 


3.72 


>30yr 


2.1 10 31 


3.78 


>30yr 


6-10- 29 


HR6416,- 








3.67 




0.82 


3.70 


0.06 


0.85 


3.51 


3.41 


1 


3.16 


1.05 


0.19 


HR 6998 


64 


14.3 


4.91 


3.19 




1 


3.19 


-0.60 


0.34 


3.02 


1.25 


1 


2.61 


1.18 


0.076 


HR 7703 


57 


13.2 


3.21 


5.08 


B 


6.5- 10" 9 


2.00 


3.64 


2.3 10 23 


2.02 


>30yr 


2.4 10 18 


2.00 


5690. 


8-10~ 17 


HR 7703,- 








1.98 




1 


2.00 


-0.01 


0.94 


1.82 


1.23 


1 


1.65 


2.87 


0.47 


4> 2 Pav 


88 


16.9 


5.72 


6.75 




0.0099 


6.70 


-1.01 


0.15 


6.01 


10.2 


0.14 


5.95 


1.64 


0.47 


HR 8323 


104 


12.8 


3.12 


3.94 


A 


0.00015 


3.71 


-1.17 


0.0003 


2.89 


1370. 


7.8- 10" 


2.89 


1300. 


MO" 9 


e Ind A 


132 


16.9 


2.25 


4.49 


P? 





3.05 


2.36 


1.6 10 23 


3.02 


>30yr 


2.7- 10 19 


2.99 


6350. 


3-10- 18 


e Ind A, 








2.99 




1.8- 10 7 


3.01 


0.01 


0.93 


2.70 


35.7 


0.0049 


2.63 


35.7 


0.0013 


HR8501 


78 


14.8 


5.47 


22.5 


B 





4.52 


17.3 





4.54 


>30yr 


4 10 48 


4.47 


7650. 


4- 10- 46 


HR 8501,- 








4.46 




0.99 


4.49 


0.11 


0.69 


4.12 


64.1 


1 


3.77 


1.61 


0.053 


HR 8883 


45 


13.3 


7.23 


63.7 


A 





64.0 


-6.90 


0.45 


32.6 


7.60 


4.3 10 9 


28.2 


13.1 


4-10- 10 


HR 8883,- 








31.7 







32.1 


-0.76 


1 


19.9 


1.58 


1.6 10 5 


18.1 


1.22 


5-10- 6 



Listed are the number of binned observations Nbm, the combined time baseline T, the mean combined noise term cr (including jitter), the combined 
weighted rms of the time series, the flag S for the probable main source of the RV variations as concluded in this work in Sect. 5 (A - activity, B - 
binary/wide stellar companion, P - planet), the x 1 -probability for fitting a constant, and the false alarm probabilities (FAP) for the other tests. Also 
listed are the weighted scatter of the residuals (i'ms stopc , rms sin , rms Kcp ) and some best-fitting parameters (slope and the periods P sin and Picep)- 



of the power values for the best-fitting sinusoid and Keplerian 
orbit Q?best and /> K ep,best) via 

Prob(p > /?best) = (1 - Pbest)^ and (11) 

Prob(/? K ep > ^Kep.best) = (1 + ^J Z /?Kep.best)(l ~ ^Kep.best)^ , 

(12) 

respectively. Compared to the probability functions given by 
these authors which account for one offset, here are slight modi- 
fications in the equations (numerator in the fractional terms de- 



creased by 2) arising from the three zero points, i.e. two more 
free parameter^. 

The final false alarm probability (FAP) for the period search 
accounts for the number of independent freq uencies M with the 
simple estimate M as fT (Cumming 2004), i.e. the frequency 



19 The corresponding normalisation as 

z= (yLstan,-* 2 )/^- Wan,) & ^ 

distribution, respectively. 
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Table 5. Information about wide companions. Listed are separation p, period estimate P, references, and our estimates for minimum 
masses M m [ n derived via Eq. (0 from the slope of the linear fit to the data (Table [4]l and p. For minimum masses in brackets the 
trends were not significant. 



Star 


Companion 


Pi"] 


P[AU] 


P[yi] 


Ref. 


M mm [M lup ] 


a For 


GJ 127 B (G7V) 


4.4 


62 


314 


[BP, H, P] 


248 


a Men 


HD 43834 B (M3.5) 


3.05 


31 




[E] 


(1) 


HR 2667 


GJ 9223 B (KOV) 


20.5 


332 




[F, WD] 


(191) 


HR 4523 


GJ 442 B (M4V) 


25.4 


234 




[P] 


(3) 


GJ 570 A 


GJ 570 BC (M1.5V+M3V) 


24.7 


146 




[B] 


325 


HR 6416 


GJ 666 B (MOV) 


10.4 


92 


550 


[LH, P] 


446 


HR 7703 


GJ 783 B (M3.5) 


7.1 


43 




[PJ 


38 


6 Ind A 


GJ 845 Bab (T1+T6) 


402.3 


1459 




[S] 


28130 


HR 8501 


GJ 853 B(V< 10 mag) 


2.5-3.4 


41 




[M, WD] 


163 



further companions 



GJ570D(T, 258'.'3) 
GJ 666 C (M6.5V, 41'.' 8) 
GJ 666 D (M7V, 40/7) 



References: [BP] Baize & Petit! 1989), [B] Burgasser et al. (2000), [E] Eggenberger et al. (2007), [F] Favata et al. ( 1997), [H] Heintz ( 197 
iLuvte n & Hughes ( 1980), [M] Mason et al] ( 1200 ll) . IPl lPoveda et alj dl9941) . rSl lScholz et al.l <2003h . IWDl lWorlev & Douglass! dl99 



[LH] 



0.4 
0.3 



HR2667- 




100 1000 
Period [d] 



10000 



Figure 5. GLS (top) and Keplerian (bottom) periodogram for 
HR 2667 which is constant to 3.1 m/s. The horizontal lines mark 
the 1(T 2 , 1(T 3 , and 1(T 4 false alarm probabilities (FAP). 



range / and the time baseline T, and is given by 
FAP= l-[l-Prob(p>p best )] M 



(13) 



and can be approximated by FAP * M ■ Prob(/? > pbest) for 
FAP «c 1 . Since our frequency search interval ranges from to 
1 d _1 , we have typically M ~ 5500 for a 15 year time baseline. 

Table [4] summarises the formal best-fitting sinusoidal and 
Keplerian periods (P sm and Picep) found by the periodograms 
along with their residual weighted rms and FAPFI Our approach 
recovers all stars that exhibit long-term trends emulated by long 
periods and generally decreases the rms down to a few m/s. 

We identify for HR 506, l Hor, and HR 3259 the same 
periods that were p reviously ann ounced as planet ary signals 
ilKurster et al.1 l2000. Mayor et al@ lLovis et al.ll2006h . In Sect.|5] 
we derive for HR 506 and t Hor refined orbital solutions (see 
also Table |6] and |7]) and investigate the correlation between RV 
and activity indicators. 

For the refinement of the orbital para meters and the error 
estimation we used the program GaussFit ( Jeffe rvs et al.| [l988) 



20 We list the formal, analytic FAP for the Kepl erian orbits, but we do 
not hig hlight them in the table. As remarked in Zechmeiste r & Kiirsteil 
( 2009), this FAP is likely to be underestimated due to an underestimated 
number of independent frequencies. Also Keplerian solutions tend to fit 
outliers (likely orginating from non-Gaussian noise) making them less 
robust for period search. 



which can solve general nonlinear fit-problems by weighted least 
squares and robust estimation. As initial guess we provided the 
parameters found with the Keplerian periodogram in the previ- 
ous section. All offsets were free parameters. 

We provide also a first orbital solution for the spectroscopic 
binary HR 2400 (Table |8). However, for k For, a Cen A+B and 
the giant HR 3677 it is not possible to give a reliable orbital solu- 
tion since our measurements cover only a small piece of their 
orbits. Companion masses are estimated in Sect. 5. 

4.5. Correlations with Ca II H&K, BIS, and FWHM 

Radial velocity variations can be caused by stellar variability 
such as oscillations, granulations, spots, and magnetic activity 
cycles. They can affect the stellar line profile and, in case of 
spots and magnetic cycles, also the amount of Ca II H&K emis- 
sion. To test this, one can analyse the /?^ K -index and the shape 
of the cross-correlation function (CCF) function in particular its 
full width at half m aximum (FWHM) and bisector span (BIS, 
Oueloz et al. 2001) which are measures for the averaged stel- 
lar line width and asymmetry, respectively. The activity indic- 
ators are products of the HARPS pipeline and the c omputation 
of (Shk and) R' HK is described in lLovis et al.l d201 ll) . Note that 
these indicators do not cover the whole survey, because they can- 
not be derived from the CES spectra (which do not include the 
Ca II region and are contaminated by th e iodine l ines). The 5hk 
errors are derived from photon noise (lLovis etalj2011h . the BIS 
span errors are taken as twi cePI the internal RV errors, and the 
FWHM errors are 2.35 timeothe internal RV errors. 

The time series of the activity indicators and their corre- 
lations with RVs are shown for each star in Figs. IB.H4B.33l 
Tables IA.4I and IA.5I summarise for each star and activity indi- 
cator their mean values, scatter, and the correlation coefficients. 
Statistically significant linear correlations with FAP < 0.01 are 
highlighted in bold font in the Tables and depicted by a solid line 
in the Figures. 



21 The precision of measuring the bisector velocity in the upper and 
lower part of the CCF (i.e. each uses only a half of the gradients in the 
CCF) is ss y/2cr RV and when taking their difference adding both errors 
in quadrature yields another factor of y2, i.e. a factor of two in total for 
the BIS span. 

22 For a Gaussian fit the mean parameter errors for center (RV) an d 
width (s) are the same cr RV = cr s (e.g. Eq. (5.8) in lKaper et al.l l 1 966). 
Moreover, since for a Gaussian function FWHM = 2 V21n2- s, we have 

C"FWHM = 235(T RV . 
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Note that a high statistical correlation does not necessarily 
mean a physical correlation, in particular when both quantit- 
ies exhibit just trends which could coincide just by chance and 
temporarily. However, if the correlation is present during more 
complex temporal variations, e.g., both quantities have the same 
period, a planetary hypothesis should be excluded. But note also, 
that the Sun hosts a Jupiter in a 12 yr orbit and shows a compar- 
ably long magnetic cycle (1 1 yr). 

In the sample, t Cet has the smallest variations in logR' nK 
(0.005 dex), while f 1 Ret has the largest variations (0.048 dex). 
We find significant correlations between RV and log R' HK for the 
stars yS Hyi, HR 209, ( 1 Ret, a Men, HR 4979, HR 8323, and for 
the two stars with planet candidates l Hor and e Ind A. However, 
the correlation seen for a Cen B (probably also HR 7703 and 
e Ind A) is artificial since the RV trend is largely caused by 
a wide companion, respectively, instead of a magnetic cycle 
(Se ct. 5). 

lLovis et ail (1201 ll) provide also a relation to estimate the 
slope of the RV and R' nK correlations based on the stellar tem- 
perature T e g and metallicity [Fe/H]. After conversior0 to a slope 
w.r.t. logR' UK by multiplication with a factor of In 10 • R' HK /10 5 , 
these estimates can be compared to the derived slopes given in 
Table IA.4I As an additional cross-check, we will do this com- 
parison occasionally in Sect. [5] when we conclude for a magnetic 
cycle hypothesis. We no te that eight stars^ n were included also 
in a sample analysed by lLovis et al.l (1201 ll) for magnetic cycles 
via With the exception of the RV standard star t Cet and 
the subgiant 6 Eri, these authors reported R' HK cycles/trends for 
these stars. 

Significant correlations between RV and BIS are found for 
7 stars (t Cet, t Hor, e Eri, a Men, HR 4979, HR 8323, and 
HR 8883). The correlation for three more stars (a Cen B, 
GJ 570 A, and HR 8501) should be artificial, since the RV trends 
can be attributed to a wide stellar companion. 

Finally, RV-FWHM correlations are found for 7 stars (J3 Hyi, 
HR 209, t Cet, ? Ret, Ret, a Men, and HR 8323). For 6 
other stars (HR 2400, HR 3677, a Cen A, a Cen B, GJ 570 A, 
and probably also e Ind A) the correlations are artificial due the 
RV trends caused by their wide companions. 

RV variations caused by the magnetic cy cle should result i n 
positive correlations with all three indicators (lLovis et alj|201 lb . 
The stars a Men and HR 8323 are nice showcase examples for 
this effect. On the other hand, an RV-BIS anti-correlation (cf. the 
active stars i Hor, e Eri, and also the g iant HR 8883) is expected 
for rotating spots (Bois se et alJl201 11) and should be therefore 
related to the stellar rotation period. 

4.6. Detection limits 

To demonstrate the sensitivity of our survey, we have calcu- 
lated for each star conservative 99.9% detec tion limits for cir- 
cular orbits following the method outlined in Zechm eister et al.l 
(2009). As an example, Fig.|6]illustrates the upper mass limit for 
HR 2667 (one of our most constant stars) showing that we are 
sensitive approximately to Jupiter analogues. Because the more 
precise HARPS data typically cover only 1500 d, there is a loss 
of sensitivity for longer periods indicated by a steep increase of 



Orbital Period P [d] 
10 100 1000 10 4 



In 10-; 



: In 10 • RHK-fj^ = In 10 • RHK ■ C RV 
C RV (T ef f, [Fe/H]) is given by Eq. (9) in 



AlagRHK AinRHK 

where the sensitiv ity Crv 
lLovis etal.U20Uh . 
24 f Tuc (HD1581), r Cet (HD 10700), f Ret (HD 20807), 6 Eri 
(HD 23249), HR 3259 (HD 69830), HR 4979 (HD 114613), HR 8323 
(HD 207129), and e Ind A (HD 209100) 




0.01 



0.01 



0.1 1 10 

Orbital Distance a [AU] 



Figure 6. Detection limit for HR 2667 considering circular or- 
bits. The cross marks the distance and m sin i of a Jupiter ana- 
logue for i = 90°. The vertical lines indicate the time baseline of 
the HARPS and all combined measurements, respectively. 



1000 




0.01 



0.1 1 

Orbital Distance a [AU] 



Figure 7. Sensitivity of the survey. For each orbital distance the 
minimum mass is indicated that could have been detected for 4, 
19, and 28 stars of the 31 stars (corresponding to a fraction of 
13%, 61%, and 90%) with 99.9% significance. The cross marks 
again a Jupiter analogue. 



the upper mass limit. The longer time baseline gained with the 
CES data pushes a bit down the limit at longer periods. 

The detection limits of the other stars have a qualitatively 
similar shape to that shown in Fig. [6] For four stars (£ Tuc, t Cet, 
and the residuals of HR 7703 as well as e Ind A) the upper mass 
limit is lower than 1 Mj up at 5 AU (due to the lower stellar mass 
of 0.7 M ). For 19 stars the limit is still lower than 2Mj up and 
for 28 stars lower than 4Mj up at 5 AU (see Figure[7]). 

5. Discussion on individual targets 

In this section we discuss individually those stars that exhibit 
variability. 



P Hyi: For p Hyi lEndl et al.l d2002l) announced a trend of 
7m/s/yr with a remaining scatter of 19m/s. Here the best 
common trend is only 1 .79 m/s/yr depicted with a black solid 
line in Fig.[T9l(plus the secular acceleration of 0.86 m/s). The 
scatter around the HARPS data decreases to only 2.3 m/s (see 
Table IA.U . However, the trend increases the VLC scatter 
from 7.4 m/s to 9.0 m/s and the fitted LC-VLC offset departs 
by 2.2cr (-17 m/s, Table lA.3"l l from the measured offset. Thus, 
it is unclear whether the trend is steady. 
Additionally, we plot the 4300 d period tabulated for p Hyi 
(Table [4]i with a black dashed line in Fig. [19] This period 
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Table 6. Orbital parameters for the planetary companion to 
HR 506. 



Parameter 




Value 


P 


[d] 


994.2 + 8.6 


K 


[m/s] 


17.3 ± 1.0 


To 


[JD] 


2 450 088 ±25 


to 


[°] 


(fixed) 


e 




(fixed) 


a 


[AU] 


2.05 + 0.24 


M sin i 


[M Jup ] 


0.94 ± 0.05 


N 




158 


rms 


[m/s] 


7.8 



matches that of a Jupiter analogue, while the amplitude of 
6.5 m/s would result in a formal minimum mass of 0.56 Mj up . 
Compared to the trend in the previous section the fitted LC- 
VLC offset is less discrepant (-8.0 m/s), but the sine fit is 
less significant than trend and still not supported by the 
VLC data, because their scatter increases from 7.4 m/s to 
8.9 m/s. Moreover, logR' HK and FWHM correlate with the 
RVs. Hence the long-term variations might be related to the 
magnetic cycle. 

HR 209: logR' HK and FWHM correlate with the RVs. Hence the 
RV variations are related to stellar activity probably induced 
by spots. 

HR 506 (HD 10647): A planet candidate was presented 
by Mayor et aP based on CORALIE measurements. 
iJones et alJ ((2004) found also weak evidence for a similar 
signal with t he AAT, but did n ot exclude stellar activity 
as the cause. iButler et all d2006l) listed AAT RV data and 
derived orbital parameters. 

For HR 506 we clearly recovered the long RV period in the 
period analysis. Hence, we combine our observations with 
AAT data and CORALIE data to fit the orbit. Because more 
cycles have been covered, our combined solution gives a 
more precise period compared to the solutions given by the 
other authors (P = (1003 + 56)d, e = 0.16 + 0.22 and 
P = (1040 ± 37) d, e = 0.18 ± 0.08, respectively). For 
our three combined data sets an eccentric orbit does not 
fit much better (Table |4| and also in the solution for the 
five combined data sets the eccentricity vanished. Therefore 
a circular obit was fitted (e and a> fixed to zero, Table [6j 
Fig. [8]). The semi-axis and the companion minimum-mass 
were derived by assuming a stellar mass of 1.17 ± 0.1 M m 
(TableQ}. 

There is no clear RV-log/?^ correlation (r = 0.45, 
FAP = 1.7%). However, when subtracting the 994 d 
RV period, there is a significant correlation (r = 0.63, 
FAP = 0.033%), implying that the residuals of this act- 
ive star are affected by stellar activity. The RV-BIS and 
RV-FWHM correlations are not significant, also not for 
the residual RVs. The RV residuals of HR 506 also have 
some excess variability and a marginally significant trend 
of -2.41 m/s/yr which however increases the rms of the LC 
data (from 18.8m/s to 20.0m/s, Table[0}- 
t Cet: The small trend reduces the scatter from 1.37 m/s to just 
1.33 m/s and is just significant due to the large number of 
observations. However, at this level an instrumental cause is 
likely for the small trend and the RV-BIS and RV-FWHM 
correlation. In this respect we also like to point out that the 
FWHM of r Cet (but also some others stars, e.g. £ Tuc) ex- 
hibits a noticeable positive long-term trend (Fig. lB.7| ) which 
might be due to a drifting focus of HARPS. In this case, as- 



suming a constant line equivalent width, a negative trend is 
expected and indeed see n in the contrast (depth) of t he CCF 
(recently also noted by iGomes da Silva et alj|2012l) . t Cet 
has the smallest variations in log/? HK in the sample. 
k For: The RVs of k For decline over the whole time baseline of 
14 yr which indicate s an orbital peri od longer than the estim- 
ate of 21 yr given in End l et all (120021) . The Keplerian period 
of 10700 d listed in Table|4](29.3 yr) is not well constrained. 
However, again with the slope and Eq. (0 this period might 
be used to assess a minimum mass of 0.36 M for the com- 
panion. The RV residuals do not exhibit significant variabil- 
ity. 

l Hor: For this active star IKurster et al.l (120001) discov ered a 
planet. The signal was also seen by iNaef et al.l (12 001 ) using 
the CORALIE spectrograph and bv lButler et alj d2001l) with 
the AAT. Using the HARPS data of IVauclair et alJ d2 008) 
taken for an asteroseismology campaign, iBoisse et al.ll201 ll 
searched for short period companions (P < 7 d), but detected 
no further planet. 

In our analy sis we recovered th e planetary signal seen 
for l Hor by IKurster et ail d2000t) . We refined the orbital 
parameters with ou r data sets and inclu ded also the AAT 
(IButler et all 120061) and CORALIE data dNaef et all l200ll) 
as shown in Fig. [9] For these two data sets the jitter es- 
timate (from Table IA.2b was added to the measurements 
error. Due to the activity of the star the residuals have a 
high scatter. The orbital parameters are listed in Table [7] 
The semi-axis and the companion minimum-mass were de- 
rived by assuming a stellar mass of 1 .25 ± 0.02 M m (TableQ}. 
An astrometric upper mass lim it of 18.4Mj up was placed by 
lReffert& Ouirrenbachl d201ll) in c ombination with the or- 
bital solution of Butler et al. (2006). The orbital period of 
307 d is different from the rotation period of ~8 d (see be- 
low), while a relation with the magnetic cycle has to be dis- 
cussed 

iMetoalfe et al.l d2010l) reported that a magnetic activity cycle 
of ~ 1 .6 yr bears no obvious relation to the orbital period, but 
we find that there is a moderately significant RV-log/?^ K 
correlation (Fig. |B.9t and that the refined orbital period 
(307 d) is half of the formal best-fitting period (630 d) for 
the combined HARPS and SMARTS S HK-index measure- 
ments (Fig . [TOl we us e the S HK-index instea d of logR' nK for 
the comparison, since Metcalfe et al. (2010) published 5 HK- 
index measurements). However, if the 307 d RV period were 
be caused by a magnetic cycle, we would expect the same 
period f or the S HK-index and a positive correlation with 
the BIS dLovis et alj|201 ll) , while there is a formal RV-BIS 
anti-correlation probably artificially induced by the short- 
term variations (Fig. ITTT) . Therefore we tend to prefer the 
planet hypothesis as a more simple and likely explanation. 
Moreover, after subtracting the 307 d period, the correlation 
with log/? HK remains and the anti-correlation with BIS in- 
creases and becomes more significant. 
The periodogram of the residuals for l Hor (Fig. IT2l shows 
power at periods of 5.7 d as well as 7.94 d and 8.45 d (also 
found bv lBoisse etal]|2011l based on HARPS RVs). The 
latter peak coincide s with periodic variati ons (8.5 d) found 
in the Shk index by [Metcalfe et alJ d2010h . which therefore 
probably indicates the rotation of this star. The residual RVs 
anti-c orrelate with BIS w hich is indeed expected for rotating 
spots dBoisse et alJl201lh . 
a For: GJ 127 B can explain the trend of -11.5 m/s (Sect. 1431 
Table|5] Fig. [24]). So far a trend is a sufficient model. The RV 
residuals do not exhibit significant variability. 
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Figure 8. Left: RV time series for HR 506 combined with AAT and CORALIE data. Right: RVs phase folded to the period of 
P = 995 d and the residuals (bottom). 



Table 7. 

i Hor. 



Orbital parameters for the planetary companion to 



Parameter 




Value 




P 


[d] 


307.2 ± 


0.3 


K 


[m/s] 


65.3 ± 


2.2 


To 


[JD] 


2449 110 ± 


9 


CO 


[°] 


35 ± 


10 


e 




0.18± 


0.03 


a 


[AU] 


0.96 ± 


0.05 


M sin i 


[M Jup ] 


2.48 ± 


0.08 


N 




205 




rms 


[m/s] 


14.5 





( l Ret: The relatively large scatter of 12.5 m/s is indicative of 
excess variability and the periodogram of the RVs shows 
peaks close to our FAP threshold at 21.5 d and 426 d (Fig.[T"3l> 
with the first giving the lowest for the combined rms value 
(but increasing the rms of the VLC data). The strong RV- 
log R' HK (also RV-FWHM) correlation shows that most of the 
scatter is due to stellar activity, likely not only the magnetic 
cycle, since the observed correlation slope 219 + 24m/s/dex 
exceeds the predicted valuel40 ± 7m/s/dex the for a mag- 
netic cycle (Sect. |4"3] assuming 5746 + 2 7 K and [Fe/H] = 
-0.22 ± 0.05 from del Peloso et"aTJl2000h . r = 0.88 indic- 
ates that th e subtr action of this correlation will reduce the 
scatter by Vl - r 2 to 5.9 m/s. ( 1 Ret has the largest log /?'- 
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4000 4500 
JD - 2450000 



5000 



5500 



Figure 10. S-index measurements for t Hor. iMetcalfe et al.l 
(2010) published the SMARTS data (blue open circles) and the 
566 d period (dotted line). The formal best fit for the combined 
data set is 630 d (dashed line). 



variations and a quite high jitter estimate of 3.8 m/s within 
our sample. 

£ 2 Ret: Given the very simila r stellar paramet e rs o f the bin- 
ary pair Ret and £ 2 Ret, Idel Peloso et al.l J2000) already 
pointed out the baffling fact of their very different activity 
level. We do not find any significant RV excess variability 
for £ 2 Ret. However, there is a significant RV-FHWM correl- 
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Figure 9. (Left) RV time series for t Hor combined with AAT and CORALIE data. (Right) RVs phase folded to the period of 
P = 307 d and the residuals (bottom). 



ation and all three activity indicators seem to exhibit a small 
trend, which might be either due to a magnetic cycle or in- 
stru mental (cf. r Cet). 
eEri: lHatzes et al.l (l2000h announced a planet in an eccent- 
ric orb it (e = 0.6) around this active star. iBenedict et al.l 
(2006) refined the orbital solution and combined the RVs 
with astrometric observations with the HST Fine Guidance 
Sensor indicating an or b ital in clination of i = 30°. Likewise, 
Reffert & Ouirrenbach| (I201 ll) derived a similar value for 
the inclination and an upper mass limit of 6.1Mj up using 
Hipparcos astrometry and the RV orbit solution. 
However, in the periodogram for e Eri (Fig. fT4l we cannot 
find any e vidence for the long -period planet (P=2500d) sug- 
gested by Hatz es et al.l ((2000) whose orbital solution is plot- 
ted in Fig. [20] for comparison. Due to the higher precision 
of our data, the combined rms is 8.2 m/s for fitting a con- 
stant, while lHatzes et al.l (|2000|) list an rms ranging from 1 1 
to 22 m/s for different instruments in the residuals of their 
Keplerian model, e Eri is an active star and all three activity 
indicators are variable, and the BIS shows an anti-correlation 
with RV The stellar r otation period is 1 1 .2 d ( Donahu e et all 
Il996t lFr6hlichll2007l) . On a short time scale of 85min its 
variability is 0.86 m/s (Table fA/JJ, while the long term jitter 
estimate is 3.6 m/s (Table |A~2] |. Our best-fitting sine function 
has a period of 3.11 d, reduces the scatter from 8.2 m/s to 
only 6.9 m/s, but is not significant. We note t hat, combining 
the RVs of HARPS and all the RVs given in IBenedict et al.l 



(2006),Anglada -Escude & Butler! (120121) could also not con- 
firm the planet solution, though, without evaluating the sig- 
nificance, they suggested another best fitting, similar long- 
period; 

a Men: Eggenberger et al. (2007) reported a companion for 
a Men which is likely an M3.5 -M6.5 dwarf that was seen at 
a separation of 3.05" with the VLT+NACO. They also men- 
tioned (but did not quantify) a radial velocity drift for a Men 
seen with the CORALIE spectrograph and suggested that the 
companion causes this trend. 

In our data however, the star has no significant excess vari- 
ability (rms = 3.1 m/s) or trend. Therefore we cannot con- 
firm an RV trend for a Men, whose value should be un- 
changed over years if induced by an M dwarf with a separa- 
tio n of 3". However, since n o value for the trend was given 
by Eggenberge r et al.l (120071) . a quantitative comparison can- 
not be made. 

Moreover, the RVs correlate with all three activity indic- 
ators. The observed RV-log7?[j K slope of 57 ± lOm/s/dex 
(Table IA.4t agrees with the predicted sensitivity of Crv = 
69+4 m/s/dex (a ssuming T eS = 559 4 ± 36 K and [Fe/H] = 
0.10 ± 0.05 from ISantos e"tai1l2004l) indicating a magnetic 
cycle as cause of the RV variations. 
HR 2400: The orbital solution for HR 2400 is given in Table [8] 
Assuming a stellar mass of 1.2 ± 0.1 M e (Table [TJ the 
companion has a minimum mass of O.17M indicating an 
M dwarf at a separation of 9.6 AU (a/26). Hence, HR 2400 is 
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Figure 11. RV vs. BIS for t Hor. (Top) The linear correlation 
coefficient is r — -0.53, but the correlation appears more com- 
plex and loop-like. (Bottom) After subtracting the 307 d period 
from the RVs, r is -0.56. 
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Figure 12. GLS periodogram on the residuals for i Hor. In this 
and the subsequent figures (Figs. [T3lfT5l l the horizontal dashed 
lines depict FAP levels as in Fig. [5] 



a single-lined spectroscopic binary (SB1). Figure[24]displays 
the RV residuals for HR 2400 which do not have a signific- 
ant variability. In the FWHM there seems to be a trend which 
might be due to stellar activity or spectral contamination by 
the companion. In the residuals the RV-FWHM correlation 
vanished. 



HR 3259: For HR 3259 lLovis et al.1 (Hb06) reported a planet- 
ary system with three Neptune-mass planets with periods 
(amplitudes) of 8.667d (3.51 m/s), 31.56d (2.66 m/s), and 
197 d (2.20 m/s). All three periods are seen in the periodo- 
gram (Fig.[T31) and the corresponding peak heights are in the 
same order as suggested by the given amplitudes. We refer 
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Figure 13. GLS periodogram for f 1 Ret. 
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Figure 14. GLS (top) and Keplerian (bottom) periodogram for 
e Eri. There is no significant power at the period of the putative 
planet eErib(P = 2500 d). 

Table 8. Orbital parameters for the companion to HR 2400. 



Parameter 



Value 



P 
K 

T 

(0 

e 
a 

M sin i 

N 

rms 



[d] 

[m/s] 

[JD] 

[°] 

[AU] 
[M Q ] 

[m/s] 



9346 ±554 
1717 ±83 
2451 881 ± 16 
279 ± 1 
0.58 ±0.01 

9.6 ± 1.1 
0.17±0.01 

77 

4.8 
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Figure 15. GLS periodogram for HR 3259. The vertical ar- 
rows indicate the periods for the three planets announced by 
lLovis et al.1 d2006l) . 



to Anglada-Escude & Butler|j20T2h for a recent and detailed 
analysis of the HARPS data and updated orbital paramet- 
ers of the short period planets. We only note that with our 
own and archive (but without GTO) HARPS data alone we 
also have seen excess power at these periods, although these 
peaks did not reach our FAP threshold due to too few obser- 
vations (only 18 binned HARPS data points) and consider- 
able spectral leakage (aliasing). 
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Figure 16. RV data for HR 4523 phase folded to 122.1 d. The 
dashed curve is an orbital solution proposed by iTinnev et alJ 
d2011l) based on the UCLES data. HARPS RVs are shown for 
comparison. 



HR 3677: The orbit of the companion to the giant HR 3677 is 
eccentric indicated by the much lower residuals when fit- 
ting a Keplerian orbit (8.7 m/s) compared to a circular or- 
bit (35 m/s). Using the slope and Eq. (0 a raw estimate for 
the minimum companion mass is 0.69 M Q . The significant 
excess variability of the RV residuals (Fig. l24l is likely ex- 
plained by an underestimated jitter of this giant. The RV- 
FWH M correlation d isap pears in the residuals. 

HR4523: ITinnev et aT] d2011l) reported a Neptune-like planet 
with 16 M e minimum mass in a 122 d orbit. 
However, the periodogram of our data shows no power at 
periods of 1 22 d. Therefore we cannot confirm the planet 
reported by ITinnev et alJ d201 ll) . Their 145 measurements 
over 7yr with UCLES have a scatter of rms = 3.0 m/s. 
For comparison, the 62 HARPS RVs over 4.4 yr have also 
rms = 2.9 m/s. Figure [16] shows the RV data phase folded to 
the period of the proposed planet. When subtracting the pro- 
posed orbital solution, the scatter decreases to 2.56 m/s for 
UCLES, but increases to 3.29 m/s for HARPS. 

HR 4979: Most of the RV variations is likely caused by a mag- 
netic cycle, because the RV correlates with log R' RK and BIS. 
Moreover, the observed correlation slope is 81 + 12m/s/dex 
which is of the order of the value expected for magnetic 
cycles (65±2 m/s/dex; Sect. 14.51 assuming 7W - 5 729± 17 K 
and [Fe/H] = 0.19 ± 01 fromlSousa et al.ll2008h . 

a Cen A+B: We refer to iPourbaix et al.1 d2002l) for the orbital 
solution for a Cen A+B with a period of 80 yr. We would 
just like to point out that the fitted trends in Sect. 14.31 al- 
though obviously not a sufficient model, can be interpreted 
as a mean acceleration and that the ratio of these slopes is a 
measure of the mass ratio ft = = jjftgfe = 1.180 

which agrees with the value of ¥r — nm , ./ = 1.183 de- 
rived from lPourbaix et al.1 (120021) . 

Neither the residuals of a Cen A nor a Cen B exhibit a sig- 
nificant variability. However, using HARPS measure ments 
and a very complex analysis, iDu musaue etld] d2012l) have 
recently announced a planet candidate for a Cen B with a 
very small amplitude (0.51 m/s, 3.236 d). Due to the lower 
number of the (complementary) HARPS measurements in 
this work (21 vs. 459) we are not sensitive to this planet. 
The RV-log/^ K and RV-BIS correlations listed for a Cen B 
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Figure 17. The time behaviour of log R' UK that shows the mag- 
netic cycle of a Cen B. 

are only formally significant, since obviously the RV trend is 
caused by a Cen A and not by a stellar activity, and vanishes 
after subtracting the orbit. The trend seen in log R' nK is part 
of a longer magnetic cycle which be seen m ore clearly when 
adding recently published HARPS data by Dumusqu e et al.l 
J2012h . We estimate a period of P ~ 2420 d (Fig. [17} for the 
magnetic cycle, but a true periodi city is not secured, since 
only one cycle is covered and e.g. Bu ccino & Mauasl d2008l) 
reported a period of 3061 d with a FAP of 24%. 

GJ 570 A: The trend has a smaller FAP compared to sinus- 
oidal and Keplerian model and is therefore a sufficient model 
for the RV. GJ 570 BC can explain the trend of -2.7 m/s 
(Sect. 14.31 Table [5}. Also, there is a marginally significant 
linear correlation of the RVs with BIS and FWHM and there- 
fore stellar activity could contribute to the trend as well. 

HR 6416: GJ 666 B can explain the trend of 9.4 m/s (Sect. 1551 
Table [5}. A trend is a sufficient model. The residuals do not 
exhibit significant excess variability. 

HR 7703: A trend is a sufficient model and the value of 3.6 m/s 
(Sect. 1551 Table[5} can be explained by GJ 783 B. Therefore 
the formally significant RV-log ^^correlation should be ar- 
tificial, however an additional (positive or negative) contri- 
bution to the trend by activity cannot be excluded. After re- 
moving the trend, there is no significant correlation or excess 
variability. 

HR 8323: The listed significant period of 1372 d is caused by a 
stellar magnetic cycle. All three activity indicators show si- 
milar variations and strong correlations with RVs (Fig. IB.30l 
Table [PI and |A"3V The observed slope of 98 ± 10 m/s/dex 
(Table lA~4l i agrees well with the predicted RV-log R' HK sens- 
itivity of Crv = 108 + 2 m/s/dex for a magnetic cycle (as- 
suming 7W _ 59 37 + 13 K and [Fe/H] = 0.00 + 0.01 from 
ISousa et afll2008l) . 

e I nd A: The trend for e Ind A originally announced by 
lEndl et al.l (120021) is probably of highest interest, since it 
might be caused by a planet. The wide bin ary brown dwarf 
companion e Ind B (M Ba +M Bb = 1 12 M Jup , King et alJ2010h 
is too far from e Ind A (Table[5| Sect l4.3b . It can induce only 
an acceleration of 0.009 m/s/yr and hence cannot explain the 
measured trend. The common slope of 2.4 m/s/yr is predom- 
inantly based on the more precise HARPS data. This trend 
improves the individual rms of the LC and VLC only slightly 
(see Table IA.1I ). however the rms for both data sets is close 
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to or below the stated precision. The scatter of the VLC data 
around the trend is even onl y 4 .4 m/s, making it the most pre- 
cise star in the VLC data seQ In the overall picture the trend 
seems to be present for the whole observation. The fit places 
the offset between LC and VLC at -4.0 m/s (Table IA.3I vis- 
ible as the small jump in the fit in Fig. [23] rj = 0.52), i.e. is 
consistent with our physically estimated offset. 
With in our sample, e Ind A exhibits comparably large vari- 
ations in logR' HK (rmsRHK = 0.032 dex, Table lA"4l i. Since 
both the RVs and log R' HK (as well as FWHM) values have a 
long-term trend during the HARPS observations (Fig. IB.3lT l. 
this results in a significant correlation (r = 0.60). However, 
both trends could coincide just by chance temporarily, so that 
it is still not clear whether they are associated with the same 
long-term period (if any). Note also that the last XogR'^ ob- 
servations seem to depart from a steady trend in contrast to 
the RVs. 

HR 8501: A trend is a sufficient model and GJ 853 B can ex- 
plain its value of 17.3 m/s (Sect. 03] Table 0. The RV-BIS 
correlation should be artificial and vanishes in the residuals 
which also do not exhibit significant variability. 

HR 8883: The period of 7.6 d found for the giant HR 8883 
improves only the HAR PS rms, but not the LC and VLC 
rms (Table |AJ)- lEndl et all (|2002|) found strong Ca II 
H&K emission in a FEROS spectrum (cf. Shk = 0.533 
in Table IA.21 > and mentioned the high X-ray luminosity. 
Furthermore, the star is indicated in the Hipparcos catalogue 
dPerrvman et alJ 1 19971) as a photometrically variable star 
(variability flag 'LP, i.e. unresolved variable, e.g. irregular or 
semi-irregular variables). Finally, the large variation of the 
bisector span and the striking anti-correlation (r = -0.96) 
with the RVs as seen in the HARPS data (Fig.[T8j reveal that 
the variations and periods of this star are intrinsic to the giant 
(e.g. rotation, activity, or solar-like oscillations^!). The RVs 
do n ot con -elate with log R' HK or FWHM (Table lA.4l and lA.5l 
Fig.|B33j. 



6. Summary and Conclusion 

Our sample consists of three planet hosting stars. These planets 
are already known and we are able to trace their periodic RV 
signals allowing confirmation of their rea lity. The planet ho sting 
stars are t Hor, discovered by this survey (iKiirster et al.l200 0). as 
well as HR 506 (Mayor et alJ3) and HR 3259 (lLovis etaiT 2006) 
which even has a system of three planets. However, we have 
no clear additional planet detection in our sample, while some 
promising or controversial cases like /3 Hyi, e Ind A, or e Eri call 

for follow-up observations. 

Compared to our previous results in lEndl et all (120021) . we 
confirm that: 

- k For, HR 2400, and HR 3677 are SB1 spectroscopic binaries 
with long periods. 

- a For, HR 6416, and HR 8501 exhibits trends consistent with 
their wide stellar companions. 

25 6 Ind A has the sp ectral type K which has numerous stellar lines 
(high quality factor O. lBouchv et alj |2001) allowing more precise RV 
measurements. 

26 With a mass estimate of 2.1 M (TablefTJ and a luminosity estimate 
of L = 60 L from visual magnitude V, distance (TablefTJ, and bolomet- 
ric correction -0.34 mag (Cox 2000), we obtain with scaling relations 
for solar-like oscillations (Kjeldsen & Bedding 1995) an amplitude of 
v osc = 0.234m/s--^-~7m/s. 
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Figure 18. Correlation between bisector span (BIS) and RV for 
the giant HR 8883. 

- e Ind A has a steady long-term trend still explained by a plan- 
etary companion. 

Moreover, we see trends in HR 7703 and GJ 570 A which 
can also be explained by the known wide stellar companions 
listed in Tabl e f51 However, th e RV trend of 7m/s/yr repor- 
ted for (3 Hyi (Endl et al] |2002l) is only 1.8m/s/yr in the com- 
bined data set and is not seen i n the VLC data. We c annot 
confirm a linear drift reported by Eggenbergeret al. (2007) for 
a Men, nei ther can we confir m the planet HR 4523 b an- 
nounc ed bv iTinnev et alJ (1201 ll) . nor the planet e Eri b repor- 
ted bv | Hatzes et al.l (l2000l). For e Eri b, the astro metric results 
dBenedict et alj2006HReffert & Ouirrenbachl201 ll) seem to con- 
firm the planet, but they do not constitute independent detec- 
tions, since they rely on the combination with less precise RV 
measurements. Also the astrometric phase coverage is not com- 
plete. However, the same is true for our HARPS measurements. 
Hence, this planet is not yet fully disprove d with our RV non- 
detec tion or by imaging non-detections (e.g. lJanson et al.1 2007, 
2008). Furthermore, there are also discussions whether an ob- 
served inner disk of warm dust can coexist with this inner planet 
dReidemeister et al.ll2f3Tlh . 

It has to be pointed out that stellar magnetic cycles (the solar 
cycle length is 2x11 yr) may cause long-term RV variations 
which might be mistaken as companions. T he Call H&K lines 
are an appropriate indicator to check thi s dSantos et al.l 12010: 
Dumusau e et al.ll201 lat lLovis et aTll201ll) . Indeed, some stars, 
like e.g. HR 209, C Ret, and HR 8323 exhibit clear RV-log R' HK 
correlations which originates from stellar activity and prevented 
us from claiming planets from the RV variations/periods. Hence, 
also in the search for Jupiter analouges, it is important to mon- 
itor activity indicators. Unfortunately, only the 4-5 yr of HARPS 
data contain the Call H&K lines, so that the whole time baseline 
is not covered, which would have been especially desirable for 
HR 506 and t Hor. 

We have not detected a Jupiter analogue, while our upper 
mass limits demonstrate that we are sensitive to Jupiter-mass 
planets up to 5 AU, i.e. Jupiter analogues (Fig.|7]l- In our sample 
the planet candidate with the longest period is HR 506 b (2.7 yr). 
While it has a minimum mass of one Jupiter mass, its period is 
still 4 times shorter than that of Jupiter. Although our sample 
size is too small to provide a meaningful number for the occur- 
rence rate of Jupiters and Jupiter analogues, the two Jupiter-mass 
planets in our sample imply that our results are in agreement 
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with the planet frequency of othe r much larger surv eys. For ex- 
ample, from the ELODIE survey iNaef et al.l (l2005h estimated a 
fraction of 7.5 ± 1.5 % fo r stars hosting giant pla nets with peri- 
ods less than 10 yr, while ICumming et al.l d2008) der ived a fre- 
quency of 12 ± 1.6 % from the Keck survey and iMavor et al.l 
$20Ul) 9.7 ±13% (for > 100 M e ) from the HARPS/CORALIE 
survey. For Jupiter an a logs, i.e. only in the range of 3-6 AU, 
iLineweaver & Grefherf {2003) estimated a frequency of 5+2% 
which is e.g. co nsistent with the discove ry rate in the AAPS sur- 
vey (3.3+1.4 %. IWittenmver et all201 ll) . Such low observed fre- 
quencies are compatible with the non-detection of a Jupiter ana- 
logue in our sample and a decreasing frequency of giant planets 
at large r distances (> 5 AU) is also predicted with core accretion 
theory dMordasini et alj|2012l) . 

The main reason for the limited VLC precision of ~9m/s 
is the short wavelength coverage of only 39 A. Useful iodine 
lines cover a total wavelength range of 1000 A. If this 25 time 
larger range could have been used, the extrapolated VLC pre- 
cision would be 5 times higher, i.e. 1.8 m/s. The actively sta- 
bilised HARPS spectrograph has a larger wavelength coverage 
and an outperforming precision. The higher resolution of the 
CES + VLC is not necessarily an advantage because the stel- 
lar lines are nearly resolved already at the lower resolution 
of HARPS, i.e. they become not sharper and therefore do not 
provide much more RV informatioto Nevertheless the CES 
data are valuable for extending the time baseline to 15 years. 

The problem of combining long-term precision RV data from 
different instruments complicates the analysis. We could derive 
the zero point offset between LC and VLC only with a limited 
precision which leads to a loss of sensitivity for trends and long 
periods. However, this problem can occur in long-term surveys 
quite frequently as spectrographs receive upgrades or survey 
projects are transferred to new instruments. Long-term access 
to the same instrument/instrument configuration is therefore im- 
portant. 

While the RV method probes the inner region of the stel- 
lar environment, the outer regions can be explored for planets 
with direct imaging. Since both methods complement each other, 
their combination leads to a more complete picture about the ex- 
istence and nature of planets around stars. An example is e Ind A, 
where despite the high imaging sensiti vity with HST/NIC MOS 
and VLT/NACO dGeifiler et al.l 120071: Uanson et all 120091) the 
non-detection puts constraints to the companion which induces 
the observed RV trend of 2.4 m/s/yr. Hence the results and detec- 
tion limits from our survey can be valuable for other campaigns 
which target these bright stars. 

Some of our nearby and bright stars are also the subject 
of projects searching and studying surrounding debris disks. 
Structures in these disks permit conclusions about the presence 

27 The thermal line broaden ing of a line with frequency / is given by 
(e.g. lUnsold & Basche1dl2002l) 

a*™. = J^K i = _l_ HK 

f \ mc 2 1.4 ■ 10 6 ^m/amu 

and depends on the temperature T and the atomic/molecular mass m. 
Therefore resolving hydrogen lines (ra H = 1 amu) in a solar-like star 
(T = 6000 K) requires a resolution of R = {- A ~ 18 000 and for 
iron lines (m Fc = 56 amu) R = 135 000. A higher resolution as 
provided by the CES leads to an oversampling or stretching of the lines 
but does not resolve new lines or sharper features which would im- 
prove the RV measurements. Contrariwise, for iodine gas (T = 323 K, 
m\ 2 = 154 amu) the thermal line broadening effect corresponds to 
# = 2400 000. 



of outer planets. For example, the He rschel satellite ha s re- 
solved exo-Kuiper belts around <T 2 Ret dEiroa et alJ 1201 Ol) and 
around HR 506 (q 1 Eri. iLiseau et al.ll2010b . In the case of the 
RV planet-hosting star HR 506, the structure of a ring at 35- 
40 AU provides a hint for another planet. Similarly, the struc- 
ture of the known debris di sk of e Eri suggests an outer planet 
dOuillen & Thorndikell2002l) . But then for HR 8501 that has an 
RV tre nd likely due to its wide stellar companion. lEiroa et al.l 
(2010) exclude a cold debris disk. 
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Figure 19. Radial velocity time series (unbinned data). The error bars depict the internal measurement errors ARVj, i.e. not including 
jitter. LC (green open circles) and VLC (blue open diamonds) data are displayed with their measured offsets. Jumps in the curves 
indicates the difference between the measured (Sect. [331 ) and fitted (Sect.UJi LC-VLC offset. HARPS data are in red filled circles. 
The solid lines indicate significant models, while dashed lines illustrate less or non-significant alternative models. All models 
include secular acceleration. Model curves are shown for /3 Hyi (trend and long-period sinusoid, see text for discussion), and 
HR 506 (sinusoid). 
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Figure 20. Radial velocities. Continuation of Fig. [191 Model curves are shown for k For (Keplerian), i Hor (Keplerian), and a For 
(trend). The Keplerian orbit for e Eri taken from Hatzes et al] (l2000l) is not significant in this work. 
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Figure 21. Radial velocities. Continuation of Figs. [T9landl20l Model curves are shown for HR 2400 (Keplerian) and HR 3677 
(Keplerian). 
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Figure 22 . Radial v elocities. Continuation of Figs. 11914211 Model curves are shown for a Cen A and a Cen B (Keplerian orbit from 
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Figure 23. Radial velocities. Continuation of Figs. ll9l422l Model curves are shown forHR 7703 (trend), HR 8323 (sinusoid), elnd A 
(constant, i.e. only secular acceleration of 1.86 m/s/yr, and trend), and HR 8501 (trend). 



23 



Zechmeister et al.: The planet search programme at the ESO CES and HARPS 



Year 

1992 1994 1996 1998 2000 2002 2004 2006 2008 2010 



100 
80 
60 

„ 40 



o 20 



-20 
-40 
-60 



k For (GOV) 
«h CES+LC 
^ CES+VLC 
m HARPS 



49000 50000 51000 52000 53000 54000 55000 

BJD - 2 400 000 
Year 

1992 1994 1996 1998 2000 2002 2004 2006 2008 2010 
80 

60 

40 

20 



-20 

-40 

-60 

-80 

-100 

49000 50000 51000 52000 53000 54000 55000 

BJD - 2 400 000 
Year 

1992 1994 1996 1998 2000 2002 2004 2006 2008 2010 
70 



HR 3677 (G0III) 
-oh CES+LC 
CES+VLC 
HARPS 



60 
50 
40 
uT 30 

§•> 

6 10 


-10 

-20 
-30 



aCenB(KIV) i 
hh CES+LC 
h^h CES+VLC " 
i-n HARPS 




49000 50000 51000 52000 53000 54000 55000 
BJD - 2 400 000 



Year 

1992 1994 1996 1998 2000 2002 2004 2006 2008 2010 
120 




HR2400 (F8V) 
' ® ' CES+LC 
«m CES+VLC 
h^h HARPS 



100 - 
80 
60 
40 - 
20 - 
- 
-20 - 
-40 - 
-60 - 
-80 L 

49000 50000 51000 52000 53000 54000 55000 

BJD - 2 400 000 
Year 

1992 1994 1996 1998 2000 2002 2004 2006 2008 2010 

70 , 

a Cen A (G2V) 
' ^ ' CES+LC 
h^ CES+VLC 
^ HARPS 




49000 50000 51000 52000 53000 54000 55000 
BJD - 2 400 000 



Figure 24. Residual RVs (binned as described in Sect.|4|i for the spectroscopic binaries. 
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Table A.l. Residual rms detailed for each instrument resulting from the joint analysis for fitting a constant, slope, sinusoid, and 
Keplerian orbit. 
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Table A.2. Jitter estimation. <Xj it T is the weighted-^ scatter of the HARPS data in the 2h bins and r is the weightecP21 mean of 
the time coverage in these bins (which may not cover the real jitter time scales in all cases). logR' HK and activity index Shk are 
median values from the H ARPS data of this work. The jitte r estimation on long time scales £7j it long was derived as described in 
Isaacson & Fis cherl (120 1 Oh u sing B-V dPerrvman et al.lll997[) and Shr (thi s work). For comparison some literature Shk values are 
given (THl lHenrv et al.lfl996l FG1 iGrav et al . 2006, [W] Wright et al. 2004) and the jitter crji t i F io ("quadrature difference of velocity 
rms minus the formal internal errors") as given bv llsaacson & Fisch er (2010) for stars monitored with Keck. 
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Table A.3. Offset differences cvlc - clc [m/s] in the fits. For a reasonable model the expected value for the offset difference is zero 
with an uncertainty of ~ 8 m/s (Sect. [33l l. 



Star 


const 


slope 


sine 


Kep 


i Tuc 


-4.98 


-5.21 


-5.18 


-1.84 


/JHyi 


-9.18 


-16.59 


-7.96 


-8.40 


HR 209 


5.87 


-0.68 


4.87 


5.53 


vPhe 


-6.52 


-1.54 


-6.29 


-4.40 


HR506 


-9.02 


-7.02 


-6.99 


-7.71 


rCet 


-0.73 


0.22 


-0.63 


0.04 


k For 


-3010.58 


252.22 


-26.30 


-13.50 


i Hor 


8.67 


9.09 


-3.99 


-0.86 


a For 


-40.37 


-5.83 


-1.54 


-3.16 


f'Ret 


-6.12 


2.80 


-7.80 


-6.67 


C' Ret 


1.02 


-2.89 


0.65 


0.47 


e Eri 


3.06 


0.34 


3.59 


1.37 


<5Eri 


6.41 


6.90 


5.87 


6.84 


a Men 


0.49 


1.80 


0.28 


1.31 


HR 2400 


1426.64 


1830.13 


529.31 


-18.71 


HR 2667 


-0.37 


1.04 


-0.28 


-1.53 


HR 3259 


8.42 


5.30 


8.18 


8.02 


HR 3677 


2312.74 


-1760.01 


-109.17 


-5.35 


HR 4523 


0.26 


0.30 


0.15 


2.05 


HR 4979 


-11.45 


-7.81 


-11.88 


-11.90 


a Cen A 


661.11 


-83.26 


-10.67 


-0.04 


a Cen B 


-860.86 


90.80 


12.11 


-6.74 


GJ 570 A 


0.01 


8.91 


5.48 


-9.99 


HR 6416 


49.10 


3.08 


-2.41 


-2.66 


HR 6998 


-6.47 


-4.22 


-5.99 


-4.36 


HR 7703 


19.57 


2.32 


0.82 


1.96 


2 Pav 


-12.30 


-9.55 


-11.92 


-12.10 


HR 8323 


-3.18 


0.86 


-2.56 


-3.35 


e Ind A 


9.17 


-3.94 


-5.37 


-4.87 


HR 8501 


53.87 


-7.32 


-4.38 


-4.34 


HR 8883 


-13.17 


-2.41 


-12.89 


-17.12 
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Table A.4. Correlations of RV with log R' RK and BIS for the HARPS data. Listed for both are the mean, the rms, the slope, the 
correlation coefficient r and its false alarm probability. Significant correlation coefficients with FAP < 0.01 are in bold font. 



Star 


<!°g*HK> 


rmsi og RH K 


slope 


'RV-logRHK 


FAP 


<BIS> 


rms BIS 


slope 


7RV-BIS 


FAP 




[dex] 


[dex] 


[m/s/dex] 






[m/s] 


[m/s] 


[m/s/(m/s)] 






<TTuc 


-4.954 


0.008 


-20± 


14 


-0.12 


0.16 


-2.3 


0.9 


-0.15+ 0.12 


-0.11 


0.2 


ySHyi 


-5.061 


0.019 


105+ 


17 


0.64 


1.5- 10~ 7 


27.5 


1.8 


0.49+ 0.22 


0.29 


0.032 


HR 209 


-4.660 


0.029 


241± 


44 


0.82 


8.6 • 10- 5 


-0.5 


5.5 


-0.08+ 0.41 


-0.05 


0.85 


vPhe 


-4.986 


0.014 


89± 


42 


0.46 


0.05 


21.3 


1.3 


0.60+ 0.47 


0.30 


0.22 


HR506 


-4.775 


0.017 


286± 


112 


0.45 


0.017 


27.9 


3.2 


1.36+ 0.63 


0.39 


0.042 


HR 506 r 






206+ 


50 


0.63 


0.00033 






-0.56+ 0.33 


-0.31 


0.1 


rCet 


-4.953 


0.005 


-24± 


15 


-0.10 


0.099 


-46.3 


0.5 


0.48+ 0.13 


0.22 


0.00018 


k For 


-5.018 


0.032 


6126± 


3609 


0.29 


0.1 


27.8 


1.8 


88.05+63.57 


0.24 


0.18 


k For r 






-15± 


10 


-0.25 


0.17 






0.36+ 0.17 


0.35 


0.048 


( Hor 


-4.670 


0.030 


475± 


119 


0.47 


0.00019 


24.6 


9.1 


-1.75+ 0.37 


-0.53 


1.6- 10~ 5 


( Hor r 






168± 


42 


0.47 


0.00019 






-0.65+ 0.13 


-0.56 


5.4 • 10-" 


a For 


-5.009 


0.007 


-164± 


350 


-0.10 


0.64 


67.1 


2.3 


-1.50+ 1.03 


-0.29 


0.16 


a For r 






-37± 


70 


-0.11 


0.6 






0.09+ 0.21 


0.09 


0.68 


C Ret 


-4.662 


0.048 


219+ 


24 


0.88 


2.2 • IQ- 9 


-7.1 


6.9 


0.59+ 0.33 


0.35 


0.084 


<T 2 Ret 


-4.892 


0.015 


66± 


24 


0.46 


0.011 


-14.0 


1.4 


0.34+ 0.29 


0.22 


0.25 


e Eri 


-4.474 


0.024 


19± 


55 


0.07 


0.74 


27.0 


6.8 


-0.55+ 0.17 


-0.54 


0.0029 


<5Eri 


-5.201 


0.008 


-1± 


52 


-0.00 


0.99 


14.4 


0.9 


-0.05+ 0.45 


-0.02 


0.92 


a Men 


-4.959 


0.026 


57± 


11 


0.61 


5.8 • 10~ 6 


-35.1 


2.5 


0.51+ 0.12 


0.53 


0.00013 


HR 2400 


-5.086 


0.017 


2197± 


3039 


0.16 


0.48 


44.6 


3.3 


-17.42+14.96 


-0.26 


0.26 


HR 2400 r 






-54± 


33 


-0.35 


0.12 






-0.19+ 0.17 


-0.25 


0.27 


HR 2667 


-4.993 


0.007 


-60± 


38 


-0.32 


0.13 


-9.0 


0.9 


0.56+ 0.32 


0.35 


0.099 


HR 3259 


-5.005 


0.010 


-3+ 


32 


-0.01 


0.93 


-7.3 


1.1 


-0.14+ 0.29 


-0.04 


0.64 


HR 3677 


-4.812 


0.007 


-39685+28001 


-0.30 


0.17 


60.6 


2.9 


-101.83+69.70 


-0.31 


0.16 


HR 3677 r 






104+ 


230 


0.10 


0.66 






-0.33+ 0.57 


-0.13 


0.57 


HR 4523 


-4.942 


0.010 


29+ 


39 


0.10 


0.46 


-33.0 


1.2 


0.59+ 0.31 


0.24 


0.059 


HR 4979 


-5.081 


0.025 


81+ 


12 


0.55 


2.8 • 10-' 


8.8 


1.7 


1.26+ 0.17 


0.60 


5.2- 10-" 


a Cen A 


-5.063 


0.007 


-3542+ 


1958 


-0.26 


0.077 


-16.9 


0.7 


-38.38+19.36 


-0.28 


0.054 


a Cen A r 






-131+ 


39 


-0.44 


0.0018 






0.53+ 0.43 


0.18 


0.22 


aCenB 


-4.969 


0.045 


3218+ 


406 


0.88 


1.9- 10~ 7 


0.6 


2.4 


58.38+ 8.25 


0.85 


9.9 • 10~ 7 


a Cen B r 






17+ 


14 


0.27 


0.24 






0.36+ 0.26 


0.30 


0.18 


GJ 570 A 


-4.812 


0.029 


79+ 


32 


0.52 


0.022 


21.4 


2.5 


1.14+ 0.32 


0.65 


0.0024 


HR6416 


-4.999 


0.008 


511+ 


217 


0.47 


0.029 


-43.5 


0.6 


-4.83+ 3.20 


-0.32 


0.15 


HR 6416 r 






-33+ 


33 


-0.22 


0.33 






0.07+ 0.46 


0.03 


0.88 


HR 6998 


-4.874 


0.020 


-26+ 


13 


-0.41 


0.062 


-42.1 


1.0 


-0.26+ 0.26 


-0.22 


0.34 


HR 7703 


-4.994 


0.016 


-166+ 


49 


-0.57 


0.0025 


1.4 


1.4 


-0.42+ 0.67 


-0.13 


0.54 


HR 7703 r 






10+ 


13 


0.16 


0.44 






0.08+ 0.15 


0.11 


0.59 


<t> 2 Pav 


-4.969 


0.020 


53+ 


39 


0.30 


0.19 


49.0 


3.2 


-0.45+ 0.24 


-0.40 


0.076 


HR 8323 


-4.900 


0.028 


98+ 


10 


0.74 


1.8 ■ 10~ 15 


-9.3 


2.1 


1.31+ 0.13 


0.76 


3.5 • 10- 16 


6 Ind A 


-4.774 


0.032 


79+ 


13 


0.60 


7.7 • 10~ 8 


20.9 


2.5 


0.10+ 0.20 


0.06 


0.62 


6 Ind A r 






19+ 


10 


0.23 


0.059 






-0.23+ 0.12 


-0.22 


0.067 


HR 8501 


-4.936 


0.011 


-764+ 


426 


-0.40 


0.09 


-11.3 


1.9 


-6.87+ 1.91 


-0.66 


0.0022 


HR 8501 r 






65+ 


49 


0.31 


0.2 






0.54+ 0.25 


0.47 


0.04 


HR 8883 


-4.313 


0.036 


605+ 


445 


0.35 


0.2 


-184.8 


130.1 


-0.45+ 0.04 


-0.96 


2.8 • 10~ 8 


HR 8883 r 






191+ 


210 


0.24 


0.38 






-0.10+ 0.05 


-0.48 


0.074 
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Table A.5. Correlations of RV with FWHM for the HARPS data. Listed are the mean, the rms, the slope, the correlation coefficient 
r and its false alarm probability. Significant correlation coefficients with FAP < 0.01 are in bold font. 



Star 


(FWHM) 


rmspwHM 


slope 


'RV-FWHM 


FAP 




[km/s] 


[km/s] 


[m/s/(km/s)] 






4"Tuc 


7.283 


0.004 


32± 


27 


0.10 


0.23 


/3Hyi 


8.032 


0.005 


266± 


82 


0.41 


0.0021 


HR209 


7.503 


0.013 


447± 


125 


0.69 


0.003 


vPhe 


8.540 


0.005 


102± 


128 


0.19 


0.44 


HR506 


9.374 


0.011 


250± 


200 


0.24 


0.22 


HR 506 r 






187+ 


99 


0.35 


0.07 


rCet 


6.279 


0.005 


-62± 


13 


-0.28 


1.9 • 10- 6 


k For 


7.844 


0.004 


-71209+26361 


-0.44 


0.011 


k For r 






-56± 


82 


-0.12 


0.5 


i Hor 


9.905 


0.020 


353± 


198 


0.23 


0.08 


i Hor, 






207± 


67 


0.38 


0.003 


a For 


8.992 


0.005 


-113± 


465 


-0.05 


0.81 


a For r 






-40± 


93 


-0.09 


0.68 


C Ret 


7.084 


0.026 


383+ 


51 


0.84 


1.1 • io- 7 


I 2 Ret 


7.026 


0.005 


315+ 


58 


0.72 


8.8 • 10~ 6 


e Eri 


6.478 


0.029 


33± 


46 


0.14 


0.48 


<5Eri 


6.489 


0.003 


138+ 


127 


0.19 


0.29 


a Men 


6.964 


0.008 


186+ 


36 


0.61 


6.1 • 10~ 6 


HR 2400 


10.288 


0.007 


-20630+ 


5269 


-0.67 


0.00093 


HR 2400 r 






-149+ 


74 


-0.42 


0.057 


HR 2667 


6.893 


0.004 


43± 


76 


0.12 


0.58 


HR 3259 


6.152 


0.004 


-55+ 


80 


-0.06 


0.49 


HR 3677 


8.617 


0.009 


66798+17643 


0.65 


0.0012 


HR 3677 r 






188+ 


177 


0.23 


0.3 


HR 4523 


6.649 


0.004 


128+ 


88 


0.19 


0.15 


HR 4979 


7.516 


0.005 


71 + 


67 


0.11 


0.29 


a Cen A 


7.465 


0.006 


-6760+ 


1860 


-0.48 


0.00071 


a Cen A r 






-69+ 


45 


-0.22 


0.13 


aCenB 


6.282 


0.016 


6323+ 


1872 


0.61 


0.0032 


a Cen B r 






41+ 


40 


0.23 


0.32 


GJ 570 A 


6.202 


0.021 


136+ 


38 


0.65 


0.0024 


HR6416 


6.413 


0.004 


968+ 


393 


0.48 


0.023 


HR 6416 r 






66+ 


60 


0.24 


0.28 


HR 6998 


6.492 


0.005 


-87+ 


57 


-0.33 


0.14 


HR 7703 


5.773 


0.006 


293+ 


140 


0.39 


0.047 


HR 7703 r 






41+ 


33 


0.25 


0.22 


<f> 2 Pav 


9.906 


0.007 


-47+ 


112 


-0.10 


0.68 


HR 8323 


7.275 


0.005 


325+ 


69 


0.47 


1.2 • 10~ 5 


elndA 


6.167 


0.016 


177+ 


25 


0.66 


8.8 • 10" 10 


e Ind A r 






46+ 


19 


0.28 


0.02 


HR 8501 


7.034 


0.005 


1969+ 


778 


0.52 


0.021 


HR 8501 r 






218+ 


85 


0.53 


0.02 


HR 8883 


17.317 


0.236 


-53+ 


70 


-0.20 


0.46 


HR 8883 r 






23+ 


32 


0.19 


0.49 
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Appendix B: Online figures: HARPS time series and 
Correlation plots for RVs, log/?u K and BIS 
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Figure B.l. Activity indictors and correlations with HARPS RVs for £ Tuc (2h binned data, secular acceleration substracted, see 
Sect ED. 
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Figure B.2. Activity indictors and correlations with HARPS RVs for (3 Hyi (2 h binned data, secular acceleration substracted, see 
Sect. ED. 
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Figure B.3. Activity indictors and correlations with HARPS RVs for HR 209 (2 h binned data, secular acceleration substracted, see 
Sect ED. 
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Figure B.4. Activity indictors and correlations with HARPS RVs for v Phe (2 h binned data, secular acceleration substracted, see 
Sect. ED. 
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Figure B.5. Activity indictors and correlations with HARPS RVs for HR 506 (2 h binned data, secular acceleration substracted, see 
Sect ED. 
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Figure B.6. Activity indictors and correlations with HARPS RVs for the HR 506 residuals (2 h binned data, secular acceleration 
substracted, see Sect. 14. 111. 
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Figure B.7. Activity indictors and correlations with HARPS RVs for t Cet (2 h binned data, secular acceleration substracted, see 
Sect ED. 
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Figure B.8. Activity indictors and correlations with HARPS RVs for k For (2 h binned data, secular acceleration substracted, see 
Sect BID. 
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Figure B.9. Activity indictors and correlations with HARPS RVs for i Hor (2 h binned data, secular acceleration substracted, see 
Sect BID. 
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Figure B.10. Activity indictors and correlations with HARPS RVs for the i Hor residuals (2h binned data, secular acceleration 
substracted, see Sect.l4.lt. 
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Figure B.ll. Activity indictors and correlations with HARPS RVs for a For (2h binned data, secular acceleration substracted, see 
Sect.gjl 
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Figure B.12. Activity indictors and correlations with HARPS RVs for £ Ret (2 h binned data, secular acceleration substracted, see 
Sect ED. 
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Figure B.13. Activity indictors and correlations with HARPS RVs for £ 2 Ret (2 h binned data, secular acceleration substracted, see 
Sect ED. 
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Figure B.14. Activity indictors and correlations with HARPS RVs for e Eri (2h binned data, secular acceleration substracted, see 
Sect ED. 
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Figure B.15. Activity indictors and correlations with HARPS RVs for 6 Eri (2h binned data, secular acceleration substracted, see 
Sect BID. 
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Figure B.18. Activity indictors and correlations with HARPS RVs for HR 2667 (2h binned data, secular acceleration substracted, 
see Sect. EH)- 
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Figure B.19. Activity indictors and correlations with HARPS RVs for HR 3259 (2h binned data, secular acceleration substracted, 
see Sect. EH)- 
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Figure B.21. Activity indictors and correlations with HARPS RVs for HR 4523 (2h binned data, secular acceleration substracted, 
see Sect.HH. 



Zechmeister et al.: The planet search programme at the ESO CES and HARPS , Online Material p 23 



15 
10 

oT 5 



£ o 



-10 



HR 4979 





3000 3200 3400 3600 3800 4000 4200 4400 
JD - 2450000 



-5 
-5.02 
-5.04 
-5.06 



-5.08 
5 -5.1 
-5.12 
-5.14 



-5.16 



4 h 




HR4979 - 

9 

* 




* 

$ 


' V 


( 







3000 3200 3400 3600 3800 4000 4200 4400 
JD - 2450000 





15 r 




14 - 




13 - 




12 - 


CO 


11 - 


£ 


10 - 


CO 




00 


9 - 




8 \ 




7 - 




6 - 




5 t- 



HR 4979 



t'l 



3000 3200 3400 3600 3800 4000 4200 4400 
JD - 2450000 

7.55 



7.54 



.CO 

I 7.53 



7.52 



7.51 



HR 4979 



■ * 



3000 3200 3400 3600 3800 4000 4200 4400 
JD - 2450000 



12 
10 

8 
6 

oT 4 
1 2 

£ o 



^HR 4979 
81 m/s/dex • 




-2 
-4 
-6 
-8 

-5.16 -5.14 -5.12 -5.1 -5.08 -5.06 -5.04 -5.02 -5 
log R'hk 




9 10 11 12 13 14 15 
BIS [m/s] 



12 
10 

8 
6 

oT 4 
1 2 

£ o 

-2 
-4 
-6 
-8 





>* 




HR4979 






* 






































^. f # 

* 







7.51 7.52 7.53 

FWHM [km/s] 



7.54 



7.55 



Figure B.22. Activity indictors and correlations with HARPS RVs for HR 4979 (2 h binned data, secular acceleration substracted, 
see Sect. EH). 
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Figure B.23. Activity indictors and correlations with HARPS RVs for a Cen A (2 h binned data, secular acceleration substracted, 
see Sect.EjjV 
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Figure B.24. Activity indictors and correlations with HARPS RVs for a Cen B (2 h binned data, secular acceleration substracted, 
see Sect.EjjV 
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Figure B.26. Activity indictors and correlations with HARPS RVs for HR 6416 (2h binned data, secular acceleration substracted, 
see Sect. EH)- 
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Figure B.27. Activity indictors and correlations with HARPS RVs for HR 6998 (2 h binned data, secular acceleration substracted, 
see Sect. EH)- 
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Figure B.28. Activity indictors and correlations with HARPS RVs for HR 7703 (2 h binned data, secular acceleration substracted, 
see Sect. EH)- 
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Figure B.29. Activity indictors and correlations with HARPS RVs for 
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Figure B.30. Activity 
see Sect. EH)- 
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Figure B.31. Activity indictors and correlations with HARPS RVs 
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Figure B.32. Activity indictors and correlations with HARPS RVs for HR 8501 (2h binned data, secular acceleration substracted, 
see Sect. EH)- 
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Figure B.33. Activity 
see Sect. EH)- 
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